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We review here particular aspects of the connection between Eaplacian growth problems and classical 
integrable systems. In addition, we put forth a possible relation between quantum integrable systems 
and Eaplacian growth problems. Such a connection, if confirmed, has the potential to allow for a 
theoretical prediction of the fractal properties of Eaplacian growth clusters, through the representation 
theory of conformal field theory. 
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1 Introduction 


1.1 Laplacian Growth 

The phenomenon of Laplacian growth constitutes one of the basic paradigms for the appearance of 
fractal patterns in physical systems out of equilibrium. Laplacian growth models may be encountered 
in natural settings, in engineering applications and in the lab. To introduce Laplacian growth we 
describe briefly an instance of it found in the latter settings, in laboratory experiment involving a 
Hele-Shaw cell. 

A Hele-Shaw cell (see Fig. 1) is composed of two glass plates fixed in place such that a small gap, 
which may be filled with fluid, is allowed to form between them. One of the glass plates contains a 
hole. The glass plates are horizontal. To obtain a Laplacian growth pattern between the two plates a 
high viscosity fluid is allowed to entirely fill the gap between the plates. Then, fluid of low viscosity 
is injected through the hole. The two fluids are immiscible. As the low viscosity fluid is injected, a 
bubble of the low viscosity fluid expands within the ambient high viscosity fluid, a portion of which, 
in turn, escapes the Hele-Shaw cell through the cell’s perimeter. 

The interface between the low viscosity bubble and the high viscosity fluid surrounding it is seen 
to possess a intricate fractal pattern, which was intensely studied numerically and analytically. The 
interface consists of fingers on different scales, somewhat reminiscent of a snow flake, though much 
less regular. The fractal dimension of the interface has been established numerically [1] and experi¬ 
mentally [2], and is found to be, say, 1.71 ± 0.03. In addition, the multi-fractal nature of the so-called 
harmonic measure of these clusters have been studied [3,4]. 



Oil 


Figure 1: A low viscosity fluid, say air, is pumped through a hole (black circle) in one of the glass 
plates constricting the gap of the Hele-Shaw cell, which is otherwise filled with a high viscosity 
fluid, say oil. As air is pumped, oil escapes through the sides of the cell, while a droplet with an 
intricate interface is formed. Depending on the protocol of injection, the interface may become multi- 
connected. 

The rate of injection is immaterial to the shapes that appear, and as such we may assume this rate 
to be constant. Since then the area inside the interface increases linearly, we may measure time in 
units of the area of the low viscosity droplet. We may thus henceforth talk of the phsycal time and the 
area of the low viscosity droplet interchangeably. 
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A theoretical prediction of the fractal dimension and other fractal properties of the cluster has long 
been at the focus of a fairly large body of scientific research (although interest may have waxed and 
waned). The purpose of the current article is twofold. First, we wish to review an approach to tackle 
the problem based on the relation of Hele-Shaw growth to classical integrability. Second, we wish to 
propose a tentative new direction of research based on quantum integrability. 

As a definition of the fractal dimension of the interface we shall take use relation between the 
linear size of the droplet R and the area of the low viscosity droplet t. We define fhe fracfal dimension 
fhrough fhe observed power law relation befween fhe fwo: 

R^ = t. (1.1) 

We leave open fhe quesfion of how fo define of fhe linear size of fhe droplef, R. In addifion we shall 
assume fhaf fhe relation (1.1) holds sharply, namely fhaf D defined above does nof flucfuafe. 

1.2 The Hele-Shaw Prohahility Density 

One of fhe main goals wifhin fhe fheory of fhe Hele-Shaw problem is fo find fhe funcfion giving fhe 
probabilify densify fo find Hele-Shaw inlerfaces of differenl shapes. Lef us define t{C) as fhe area 
inside fhe inferface C, and fhe probabilify disfribufion funcfion , P{C), as fhe probabilify densify fo 
find an inferface C af time t{C), given some probabilify disfribufion on fhe initial condition Po{C). Lef 
be fhe fhe inferface C evolved a time Af. Namely, if C has area t, fhen has area t + Af. The 
funcfion P is defined as follows: 


P{C) = PoiC-^^^'>) (1.2) 

We shall specify fhe measure dC, wifh respecf fo which fhe probabilify densify P{C) is defined, buf 

shall only demand fhaf if is lime Iranslalion invarianl, dC = . Indeed, Eq. (1.2) makes liffle 

sense ofherwise. 

The probabilify densify P is time Iranslalion invarianl. Indeed, from (1.2) one may easily show 

P(C^*) = P(C). (1.3) 

The function P possesses scale invariance, as well. Lef us define, wifh some abuse of nolalions, P{R) 
as follows: 

P{R) = J P{C)S{R-R{C))dC, (1.4) 

where R{C) is fhe linear size of fhe interface. For Laplacian growlh we assume Eq. (1.1) holds which, 
we shall see, leads fo: 


Indeed, 


P{R)dR = dR^. 

(1.5) 

P{R) = j P{C)6 (^R - f^/^(C)) dC = 

(1.6) 

= DR^-^ J P{C)5 {R^ - t{C)) dC = 

(1.7) 

= DR^-^, 

(1.8) 
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where the last equality follows from the faet that the measure is time translation invariant. 

We note that the time translation invarianee of the probability density means that the probability 
measure ean be thought as a probabilistie uniform superposition of the interfaee at different times. 
That the probability density is also seale invariant is a demonstration of the fraetal property of the 
interfaee, namely Eq. (1.1), such the two properties together define the problem of finding fhe proba- 
bilify densify which encodes fhe fracfal dimensions. 

A furfher imporfanf commenf is fhaf in searching for fhe probabilify densify obeying scale and 
lime Iranslalion invariance, fhe problem may be somewhaf generalized in fhe following sense. The 
droplef in fhe Laplacian problem remains singly connecfed if fhe injecfion of fhe low viscosify fluid is 
always fasl enough lo rule ouf surface fension effecls. However, if fhe protocol of injection is alfered 
fhe droplef may become mulli-connecled. Imagine, for example faking a large droplef, wilh a fully- 
developed fracfal sfrucfure over many scales. The inlerface is self similar in fhaf if fhe low viscosify 
fluid is exfracfed, allowing fhe droplef to shrink, buf we re-scale fhe resulting inlerface by a proper 
faclor, one oblains fhe original droplef. Assume now fhaf before exlracling fhe low viscosify fluid 
we wail some lime and allow surface fension lo change fhe shape. This has fhe effecl of coarsening 
fhe shape on fhe lowesl scales as well as allowing for fhe possibilily for fhe droplef to become mulli- 
connecled. Then some low viscosify fluid is quickly exfracfed. We may still expecf fhaf fhe self- 
similarily, Ihough perhaps in a weaker sense, is preserved. If seems Ihen fhaf fhe single-connecfedness 
of fhe droplef should nol be imposed as a slricl condilion. As a resulf we consider probabilify densilies 
which give non-zero measure lo mulli-connecled droplels. 

1.3 Preview of the Paper 

The facl fhaf fhe lime Iranslalion invariance properly of fhe probabilify dislribulion function, Eq. (1.3), 
coexisls wilh fhe scale invariance properly, Eq. (1.5), is equivalenl to fhe scaling law R^{C) = t{C). 
We may formulate fhe problem of determining D, fhe fracfal dimension, as finding fhe scaling law 
(1.5) for a probabilify dislribulion funclion obeying fhe time Iranslalion invariance properly, Eq. (1.3). 
The currenl paper puls forward fhe proposilion fhaf candidales for fhe probabilify dislribulion function 
may be found wilhin fhe realm of quanlum inlegrable syslems and conformal field Iheory . 

Eel us briefly recounl fhe connection belween fracfal slrucfures and conformal field Iheory. The 
conlexl in which fhe Iwo concepls may be mosl successfully linked is wilhin fhe framework of equi¬ 
librium Iwo dimensional crilical phenomena. There, scale invariance, which is a consequence of fhe 
divergence of correlalion lenglhs leads lo conformal invariance, and if is assumed fhaf, al fhe crilical 
poinl, fhe slalislical mechanics problem is described by conformal field Iheory. The vacuum is a Gibbs 
slate in which fracfal objecfs, such as clusler boundaries, exisl al all lenglh scales. A fracfal cluster 
may Ihen be pinned to a parlicular poinl by applying a scaling field al fhaf poinl [5]. The dimension 
of fhe scaling field is determined by fhe power-law dependence of fhe scaling field on fhe small scale 
cul-off. Since Ihis cul-off also determines fhe probabilify fhaf fhe fracfal can be found and pinned 
lo a small box of linear size given by fhe cul-off scale, and Ihis in torn is related lo fhe notion of 
box-counting fracfal dimension, Ihere is a relation belween fhe dimension of fhe scaling field and fhe 
fracfal dimension of fhe objecl being pinned by if. 

Non-equilibrium problems, such as Eaplacian growlh, are unlikely lo be linked in fhe same manner 
lo conformal field Iheory. Indeed, if does nol seem a reasonable assumption fhaf fhe Eaplacian growlh 
inlerface may be pinned to a poinl by applying, to some vacuum conlaining already many fracfal 
slrucfures on all scales, a scaling field of conformal field Iheory. Such a scenario would mean fhaf 
fhe fully developed large area inlerface can be picked al once from a Gibbs slale wilhoul a need lo 
apply fhe nonlinear evolulion from an inilial inlerface. In olher words, such a scenario would mean 
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that there is a representation of the interfaee as an equilibrium objeet. 

Nevertheless, a eonneetion to eonformal field theory is an appealing prospeet, due to the faet that 
eonformal field fheory is inherenfly seale invarianf, jusf as fhe fraelal dimension of fhe inferfaee is, 
and due fo fhe fael fhaf sueh a eonneefion would alow for fools fo predief fhe fraelal dimensions, sinee 
eonformal field fheory nalurally supplies a diserefe speefrum of possible sealing dimensions, fhese in 
furn possibly delermining fhe fraelal dimension. 

Here we propose lhal inslead of basing fhe eonneefion fo eonformal field fheory on fhe notion 
of seale invarianee being promofed fo eonformal invarianee, a relalion fo eonformal field fheory may 
be esfablished making use of fhe fael lhal eonformal field fheory has a hidden inlegrable slruelure, a 
slruefure whieh ean be viewed as a quanlizalion of fhe inlegrable slruelure hidden in Laplaeian growlh. 
Al presenl, we are only able fo make Ihis eonneetion wilhin some limil, where only a small pari of 
fhe inferfaee, whieh happens fo be long and narrow, is foeused upon. Wilhin fhis limil, fhe inlegrable 
slruelure of Laplaeian growlh beeomes lhal of Ihe Korleweg-de Vries equation, Ihe quanlum analogue 
of whieh is Ihe inlegrable slruelure hidden wilhin eonformal field Iheory. 

To flesh oul Ibis eonneetion, we review firsl Ihe inlegrable slruelure wilhin Laplaeian growlh [6-9] 
in seelion 2. This inlegrable slruelure lurns oul lo be Ihe dispersionless limil of a reduelion of Ihe Iwo 
dimensional Toda lalliee. Sinee we are evenlually interested in Ihe quanlum inlegrable system, we 
musl diseuss how Ihis dispersionless limil is oblained from a fully quanlum system. To Ihis end 
seelion 3 and 4 deseribe how Ihe dispersionless limil is removed lo oblain Ihe full bul redueed and 
elassieal Iwo dimensional Toda lalliee. The nexl logieal step would be lo pass lo Ihe quanlum problem. 
However, before diseussing how lo quantize Ihe problem, we pass, on a elassieal level, from Ihe Iwo 
dimensional Toda lalliee lo Ihe Korleweg-de Vries equation [10]. This is done in seelion 5, and 
ineludes foeusing on a long and narrow finger of Ihe Laplaeian growlh. In seelion 6 we diseuss Ihe 
quantization of Ihe problem [11-15], whieh leads us lo Ihe quanlum Korleweg de-Vries equations. 
This quanlum problem lies al Ihe hearl of eonformal field Iheory. This fael allows us, in seelion 7, 
lo show how a eerlain objeel, familiar wilhin Ihe framework of eonformal and inlegrable field Iheory, 
may deseribe Ihe elassieal Laplaeian growlh problem, Ihrough Ihe probabilily dislribulion funelion 
P{C). Seelion 8 eonlains eoneluding remarks. 

2 Classical Integrability in Laplaeian Growth 

Here we review some of Ihe malhemalieal slruelures underlying Ihe problem of Laplaeian growlh. 
This ineludes Ihe deseriplion Ihrough eonformal maps and elassieal inlegrabilily [6,7]. 

2.1 Description Through Conformal Maps 

Lei us eonsider Ihe Hele-Shaw selup, whieh was already deseribed in Ihe inlroduelion. The ambienl, 
high viseosily fluid is deseribed by a veloeily field v. This veloeily is Ihe average veloeily aeross Ihe 
gap belween Ihe Iwo glass plates. Due lo Ihe high viseosily of Ihe fluid, one musl balanee Ihe internal 
viseous frielion foree experieneed by Ihe fluid, wilh Ihe foree supplied by Ihe pressure, negleeling 
inertia, whieh is small beeause of Ihe small separation of Ihe plates. As Ihe foree of frielion seale wilh 
Ihe average veloeily aeross Ihe gap v , wilh Ihe magnilude of Ihe viseosily, fi, and inversely wilh Ihe 
size of Ihe gap squared, ^ one oblain Darey’s law: 

'Pv^-VP. ( 2 . 1 ) 


6 



The factor 12 may be obtained by a more detailed, but basic calculation, which takes into account the 
no-slip conditions at the interface between the fluid and the glass plates, which results in a parabolic 
velocity profile. 

The dynamics are dictated by the following conditions. The motion of the interface between the 
inviscid and viscous fluid are determined by the normal velocity at the interface, which is proportional 
to the normal derivative of the pressure at the interface: 


Vn 


62 

12;u 


dnP. 


( 2 . 2 ) 


Incompressibility, V • i; = 0, requires, that there exists a real stream function, (p\ 

v = ^zxVe, (2.3) 

27r 

where 2 is the out-of-plane unit vector and Q is the volume flux, namely the volume of fluid injected 
into the cell in unit time, inserted here for convenience, as to make (p dimensionless. 

Defining the complex function v by v = Vx + iVy, and taking d, and 9 to be a derivative with 
respect to z and z, respectively {d = Q = write the relations between the 

pressure, the stream function and the velocity as: 

/,2 _ Q _ 

V =- dP = —idO (2.4) 

6/U TT 

We may define a complex analytic potential pt (where in ipt the index denotes explicitly the time 
dependence): 

(pt{z) = 0- i^^P- (2.5) 

Indeed, relation (2.4) suggests that (pt{z) is analytic dpt = 0, Taking the complex conjugate of (2.4) 
also leads, after simple algebra, to: 

y = -^Rdpf (2.6) 

27r 

The pressure is constant inside the droplets of lower viscosity, since the fluid there has neither 
inertia (due to the narrowness of the gap) nor viscosity (by assumption). We may assume this pressure 
to be zero. This means that we assume a type of dynamics in which, even after the low viscosity 
droplet breaks up into several droplets, the same pressure is maintained across the different droplets. 
This scenario cannot be attained in a Hele-Shaw cell, since only one of the droplets will be connected 
to the source of low viscosity fluid through the hole. One may imagine instead of a glass plate using 
a plate made of a material permeable to the low viscosity fluid but impregnable to the high viscosity 
fluid. The low viscosity fluid may then be supplied through the entire plate, and the different droplets 
will maintain the same pressure, which without loss of generality may be set to zero, since only 
gradients of the pressure play any role in the dynamics. 

We now consider the analytic properties of pf First note that we assume that the fluid is syphoned 
off in a radially symmetric fashion around infinity, which suggests v{r) ~ ^ as r —)■ oo. This leads 
to the following asymptotic behavior of pt : 


Mz) — i\og{z)^ as z —)• oo, 


(2.7) 
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which leads to the fact that the imaginary part of (f + t, namely the stream function (f winds by 27r 
as one encircles infinity. In fact ipt is only well-defined if a branch cufs are infroduced in fhe com¬ 
plex plane, such fhaf from each droplef a branch cuf exfends fo infinily. This can be undersfood by 
considering fhaf fhe winding by 27r of ipt around infinily has ils origin in fhe winding of 0 around fhe 
droplels. Indeed, consider fhe inlegral over fhe velocity around fhe droplef. From simple hydrody- 
n ami cal kinemafics, fhis inlegral is equal lo volume flux of fluid enlering fhe cell, on fhe one hand, 
and from fhe definilion of fhe slream function, Eq. (2.3), if is equal lo fhe inlegral over fhe gradienl of 
6 around fhe droplels: 

V X df= ^ (j) V6 ■ dr, (2.8) 

leading lo fhe conclusion fhaf fhe winding of (p around fhe droplels is 27r: 



(j) yo = 2 tt. 


(2.9) 


The Ihe variable 0 can Ihus be used as a coordinate along Ihe droplels, whose range is [0, 27r]. 

A cenlral role will be given lo Ihe inverse function of rpt, which we shall denote by fp. 

= z (2.10) 

ft{ip) is defined such lhal for (p G [0, 27r], Ihe function ft{ip) gives Ihe complex coordinate of a poinl 
on Ihe interface of Ihe droplels, wilh given 6 and al time t. The function ft is defined wilh period i27r: 

ft{ip + 27rt) = ft (ip) , ( 2 . 11 ) 

For complex p, Ihe function ft{tp) is defined by analytic continuation. As a resull of Ihe periodicity, 
Ihe function ft is defined on a cylinder, S, where poinls displaced by 27n are identified, S = C/27rzZ. 
The funclion ft{ip) has Ihe following expansion al Im((/?) —)■ —oo : 

OO 

ft{ip) = + ''' = E (2-12) 

i=-i 


which resulls from (2.7). 

Since Ihe interface is described by Ihe map ft{ip) '■ [0, 27r] —)■ z and since ft is periodic in 27rz, 
Ihere is an arbilrariness in choosing ft{ip) given Ihe shape of Ihe interface. Indeed Ihe functions ft{ip) 
and ftiiip + c)) describe Ihe same interface, Ihe only difference being lhal Ihe reference poinl on Ihe 
interface where Ihe angle 6 is considered as 0 is moved from /t(0) lo ft{^c). Thus, when describing 
Ihe dynamics of Ihe interface as Ihe dynamics of ft, one has lo bear in mind lhal Ihe dynamics of ft 
which are equivalenl lo shifting Ihe reference poinl of p are meaningless and arbilrary. Il is common 
in Ihe lileralure lo remove Ibis arbilrariness by requiring lo be real. We shall nol follow lhal 
convention, and ralher leave Ibis arbilrariness in place. 

An useful malhemalical objecl lo sludy interfaces in Iwo dimensions is is Ihe Schwarz function, 
S{z) - Ihe an analytic continuation of a function giving z on Ihe boundary. Namely for z G C, we 
have 5 ( 2 ;) = ty wilh S{z) being analytic in a neighborhood of C. One may easily ascerlain lhal Ihe 
Schwarz function can be written as 


S{z) = ft{-ipt{z)), 


(2.13) 



where, as usual, for any funetion g{z), the funetion g{z) is defined as g(z). To see that (2.13) in¬ 
deed gives the Sehwarz funetion, eonsider that for z on the boundary, (pt{z) , is purely real and so, 
ft{—i(ptiz)) = ft (lift (z)) = z. The Sehwarz funetion, S{z), obeys the ‘unitarity’ eondition: 

5(5(z)) = z. (2.14) 

This relation is proven by eonsidering it on the interfaee, where it simply follows from the faet that 
S{z) = z on the interfaee, and then trivially analytieally eontinuing the relation away from the inter¬ 
faee. 


2.2 Richardson Moments 


As was found by Riehardson [9], during the eourse of the evolution an infinite set of parameters, 
namely the Riehardson moments remain eonstant. The k’th Riehardson moment is defined as: 


f dz Adz 
^ Jv+ ‘^'TTtkz^ ’ 


(2.15) 


where denotes the domain whieh is exterior to the droplet. To see that these are eonserved eonsider 
tk- 


tk = z 


dr X V 
0 2mkz^ 


(2.16) 


We may write z ■ dfx v = | {dzv — dzv). We may re-write the right hand side of this equation using 
(2.6) as ^{dzd(p -h dzBcp) = ^dtpt, one obtains: 


Q 

SlTT^k 


dipt 


(2.17) 


As the integrand is analytie in pt, the eontour of integration may be deformed to a large eontour 
surrounding infinity in the z-plane, where the asymptotie relation pt ^ log ^ holds. This yields: 


tk 


Q 

Snr'^k 



(2.18) 


Note also that the Sehwarz funetion eneodes the harmonie moments. Indeed, Stokes theorem 
states that (2.15) ean be written as a eontour integral over the boundary of the droplets ti = 
and sinee one the eontour S{z) = z, this ean be written as: 


One may eonstruet from knowledge of the Riehardson moments, tk, and the area of the droplets, 
t its shape, given some additional data, sueh as the number of eonneeted eomponents. The Riehard¬ 
son moments are thus reminiseent of the eonserved quantities of eompletely integrable systems, sinee 
the knowledge of whieh is enough to determine eompletely the dynamies. The results of Refrs. [6-8] 
show, however, that the aetual integrable strueture behind the zero surfaee-tension Hele-Shaw problem 
treats the Riehardson moments rather as the ‘times’ eonjugate to the eonserved quantities, or ‘Hamil¬ 
tonians’, whieh may be denoted by Hk- In addition it is appropriate to treat the eomplex eonjugate 
of the Riehardson moments (whieh are themselves eomplex quantities), namely the quantities tk, as 
additional independent times, with whieh one assoeiated a set of Hamiltonians, Hk- 
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To be able to define the Hamiltonians, Hk, Hk, we must first define a Poisson braeket. This is 
given by: 


= 1 . 


( 2 . 20 ) 


With this definition of the Poisson braeket, the evolution of ft is eneoded in the following relations: 

= ( 2 . 21 ) 

whieh, due to ft{—iip) = S{ft{i(p)), ean also be written as: 

{ftM, Siftivf))} = 1 . ( 2 . 22 ) 

That the evolution of ft is determined by (2.21) ean be shown as follows. First note that (2.21) 
ean be written as: 


Imdtftdjt = 1. (2.23) 

Sinee dtft is the veloeity of a point on the interfaee with given ip, we may write it in usual veetor 
notation as v, the tilde denoting that it may not be equal to the veloeity of the fluid. Nevertheless, its 
normal eomponent must be equal to the veloeity of the fluid so we have Vn = Vn- Furthermore, d^ft 
may be written in veetor notations as Taking into aeeount that if a = + lUy and b = bx + iby 

are the eomplex notations for a and b, then the veetor produet is given by a x 6 = zlm{ab), we may 
rewrite (2.23) as: 

^ dr 

vx — = 1. (2.24) 

dip 

This relation implies Vn-^ = 1, or Now, sinee 9 and are eonjugate harmonie fune- 

tions (Eq. (2.5)), one may use the Cauehy-Riemann relations to obtain ^ yielding the 

equation governing the evolution of the droplet, Eq. (2.2). 

The evolution with respeet to the physieal time, is given by relation (2.21). Alternatively, we ean 
write: 


dtf = 


(2.25) 


showing that the designation of the Hamiltonian that generate the time translations as up is appropriate. 
We may thus write H = up. Note however, that Eq. (2.25) is almost an empty statement, in the sense 
that the evolution of / eannot be extraeted from it, as the right hand side is identieally equal to the left 
hand side irrespeetive of the time evolution of the droplet. 

Another way to deseribe the evolution is through the Sehwarz funetion. Writing the Poisson 
braeket in (2.22) explicitly and dividing by one obtains: 


dS{ft{up)) _ 1 f dftiup) dS{ftiup))\ 

dt dftM) dt dip ) ' 


Evaluating this equation niip = ipt{z) and making use of 




dip 




^ ^ dipt{z) dftiup) , dftiup) 


(2.26) 


(2.27) 
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yields: 


idipt{z) ^ diptjz) dS{ft{up)) dS{ft{i(p)) 

dz dt dip dt 


The right hand side ean be identified with sueh that we have: 

9S{z) _ idiptjz) 

dt dz ' ^ ’ 

whieh eonstitutes a simple evolution equation for the Sehwarz funetion. 

Eq. (2.29) is useful to show that the dynamies preseribed by (2.21) is sueh that the pressure is 
equal on all interfaees. This needs to be shown beeause above we have only assumed that the pressure 
is eonstant on any eonneeted eomponent of the interfaee. Suppose then that points a and b are on the 
interfaee, perhaps on different eonneeted eomponents of whieh. Eq. (2.29) together with (2.5) allows 
us to write: 


P{a)-Pib) 


d 

= Re- I S(z)dz. 


(2.30) 


On the other hand we have: 


(■h fS{h) 

/ S{z)dz = - _ S{y)dy + zS{z)\^^, (2.31) 

J a J S{a) 

whieh is obtained by first changing the integration variable from ztoy = S{z), on the left hand side, 
this change of variables being facilitated by using unitarity, Eq. (2.14), to obtain z = S{y), then 
integrating by parts the result. Taking the real part of Eq. (2.31) and a time derivative, one obtains: 

2 Re J dtS{z)dz = Ke (^zdtS{z) — zdtS{z)^ =0, (2.32) 

which when compared to (2.30) yields: 


P{a) = P{b), (2.33) 

for any two points a, b on the interface, including points on different connected components. 

2.3 Integrability of Laplacian growth 

Integrability of Eaplacian growth is the statement that the evolution of the Hele-Shaw interface can 
be described by a Hamiltonian, H, which is in involution with an infinite set of other Hamiltonians 
(conserved quantities), with respect to some Poisson bracket. As usual, the other Hamiltonians, 
in Eaplacian growth can be used to generate flows wifh respecf fo new limes tk, in one lo one cor¬ 
respondence lo Ihe Hamilfonians, H^. These Hamiltonians in Eaplacian growth can be chosen such 
that tk are the Richardson moments of the interface. Namely, the Hamiltonians, Hk, generate defor¬ 
mations of the interface corresponding to changing the Richardson moments tk, respectively. Due to 
the fact that these moments are complex, the Hamiltonians are also complex and to each Hk we may 
associate a Hk, which generates, formally, deformations which correspond to changing tk- Of course 
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one cannot change tk without changing tk, yet, as usual, they may be treated as independent variables, 
just as Hk and Hf^ may be treated as independent conserved quantities or Hamiltonians. 

The Hamiltonians read explicitly: 

= + (2.34) 

where is a an arbitrary function of the times only, which (as shall be seen below) does not affect the 
dynamics, while the index + denotes taking the positive (polynomial) part of the Laurent expansion 
around infinity in the variable e^‘^, using the following convention. Due to Eq. (2.12) we may expand 
ft as 


ft = E 

j=-oo 

The positive part of this series is defined as follows: 



(2.35) 


(2.36) 


which is natural, except the zeroth order term, defined wifh a factor 1 /2 for lafer convenience. 

The Hamilfonian Hk is defined as follows: 

Hk{^,t) = HkWJ) = + (2-37) 

We now show fhaf fhe evolufion generated by fhe Hamilfonians, 

^E^ = {/t(z(/p),i7fc}, (2.38) 

= (2.39) 

effecl a deformalion of fhe confour corresponding to changing fhe Richardson momenfs. Namely, 

we show fhaf if (2.15) holds af any momenf in limes, Ihen if holds for all subsequenf limes, for fhe 
evolufion prescribed by (2.38), (2.39). 

Firsl nole fhaf fhe complex conjugate of Eq. (2.39) can be also casl as fhe evolufion of fhe Schwarz 
function. Indeed, fhe complex equalion of fhaf equafion, has fhe following form: 

= (2.40) 

atk 

which is fhe evolufion of fhe Schwarz function, S{z) = ft{—upt{z)), evaluated al z = We 

proceed by firsl showing: 


dS{z) dHk dS{z) dHk 

dtk dz ’ dtk dz 


We do Ihis by wriling: 


dS{z) 

dtk 


dft{-iv) 

dtk 


diptjz) dft{-i(p) 
dtk dip 




(2.41) 


(2.42) 
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Using (2.38) and the first equation in (2.27) along with 


-1 


d(pt{z) 

dtk 


dftM f dft{vf) \ 
dtk \ dip ) 


‘P=Vt{z) 


one may re-write (2.42) as: 

dSjz) ^ dih ( dS{ft{iip)) dS{z) dftiiip) \ 
dtk dip V dt dz dt ) 




dS{z) dHk 
dt dip 




(2.43) 


Making use now of (2.29), the first equation in (2.41) follows. The seeond equation in (2.41) may be 
derived by the analogous manipulations. 

Given (2.41), and the definition of Hk, we may write 


dSjz) 

dtk 


kz’^-^ + 0{l/z‘^), 


dS{z) 

dtk 


0{\/z\ 


(2.44) 


Indeed, substituting into the first equation in (2.41) the explieit expression for Hk and evaluating at 
ip = ipt{z), leads to the first equation. For the seeond equation note that Hk eontains only negative 
powers of These, in turn, when evaluated ni ip = ipt{z) ean be written as negative powers of z, 
due to (2.12), leading to the seeond equation in (2.44). Equation (2.44) immediately implies: 


_idS{z) dz _ 
dtk 2mk 


^-i dSjz) dz ^ ^ 
dtk 27r^/c 


(2.45) 


Taking a derivative with respeet to tk and ik of Eq. (2.19), one realizes that Eq. (2.45) is just the 
differential form of these equations. Thus we have proven the required statement. 

Note that in (2.34) are indeed arbitrary, in the sense that we did not need to assume anything 
about them to prove that (2.19) follow from (2.38), (2.39). In faet, just generates transforma¬ 
tions that shift the referenee point for the angle ip. To see this eonsider that one always applies the 
Hamiltonians in a real eombination. Namely any deformation of the droplet is effeeted by: 


dft = '^{ft,dtkHk - StkHk]. 

k 


(2.46) 


The free ((/^-independent) term in 5tkHk — StkHk reads 2ilm{5tkC^^^), while the latter generates the 
following deformation of fp. 


{/t,2^Im((54c(^^)} = 2lm{5tkC^^^) 


dft 

dip' 


(2.47) 


where the dot signifies a t derivafive. Thus, 2ilm{5tkC^^^) indeed is fhe generator of real, times- 
dependent ip translations, and thus do not deform the interfaee. As a result, we may ehoose 
freely. 

We make the ehoiee of sueh that it is eonsistent with the involution of the Hamiltonians. 
Consider first the indefinite integral J S{z)dz = J2k^kZ^ + flog( 2 ;) + c + 0{l/z), where c is an 
arbitrary funetion of the times, the arbitrariness of c being due to the indefiniteness of the integral. We 
may ehoose c sueh that dtc = — log(ag^^), and define the result, whieh is now definite up to a time 
and z independent additive faetor as the integral of S'(z) : 


S{z)dz = '^tkz^ + t\og{z) - [ log(a[)^^)-h 0(l/z). 
fc Jo 


(2.48) 
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The integral over 5 is a generating funetion for the Hamiltonians. Indeed due to (2.41), we have: 

= + (2.49) 

Z 

for some ^-independent eonstant, We may ehoose sueh that this eonstant is zero, and write: 




dtk 


S{z)dz 


Z 


= Hk, 


dt 


S{z)dz 


Z 


= H. 


(2.50) 


This last equation (whieh is derived from the ehoiee of states that J S{z)dz is a generat¬ 
ing funetion for the Hamiltonians and, as shall be seen below, insures that the Hamiltonians are in 
involution. Indeed, The identities r,P\, = , -sPPr = aP\, , leads to: 

’ at^dti dtidtk’ dtkdti dtidtk’ 


dm^ 

_ dm 

dHk: 

_ dHi 

dti 

Z dtk 

z' dil 

z dtk 


(2.51) 


and this yields the involution of the Hamiltonians, when written in terms of derivatives at fixed ip 
rather than at fixed z; We first write: 


dm 


dti 


dHk dcptiz) 


dm 


dti dti dip 


‘P=‘Pt(z) 


dm 

dftM , 

( dftM\ 

dHk 

dti 

dti 

K d^ ) 

dip 


‘P=‘Ptiz) 


then subtraet the equation with k and I interehanged and use (2.38) to obtain: 


dm 

dti 


dm 

dtk 


{m,m] = Q- 


(2.52) 

(2.53) 


(2.54) 


The latter being the required statement of the Hamiltonians being in involution. Of eourse one may 
easily obtain the following as well: 


dm 

dti 


dHi 

dtk 


{m,Hi} = o. 


(2.55) 


3 Algebraic Solutions to Laplacian Growth 


The interfaee may be deseribed as a eurve, namely as a relation between x and y, the Cartesian 
eoordinates on the interfaee. This may be done by finding some funetion F{x, y) sueh fhaf F{x, y) = 

0 gives fhe inferfaee. The ellipse is fhe simples! example exeepf fhe eirele. For fhe ellipse we have 

~ 2 2 

F(x,y) = ^ + ^ — 1. Using eomplex eoordinafes one may alfernafively deseribe fhe eurve by a 

’’l ^’2 

relation befween z and z, for example, fhe ellipse may be deseribed by F{z, z) = F = 0, 

where F{z, z) = \ \z\‘^ + 5 -|- z^). One may solve fhis equation for z. This 

gives fhe Sehwarz funetion, S{z), whieh satisfies F{z, S{z)) = 0. If, as is frue for fhe ellipse, fhe 
funetion F{z,z) is a polynomial in bofh argumenfs, fhen ^(z), as a funetion of z, is an algebraie 
funetion. As an example, for fhe ellipse we obfain: 


5(z) 


(rf + r 2 )z + 2 x 1 X 2 ^/z"^ — rf 

22 I 2 
'1 '2 


(3.1) 
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Namely 5 ( 2 ;) contains roots and polynomials. As an algebraic function, S{z) may be defined on a 
algebraic Riemann surface. 

Defining fhe Schwarz funcfion on fhe Riemann surface arises due fo fhe facf fhaf S{z) is generi- 
cally a mulfi-valued funcfion, such as fhe funcfion in (3.1), where fhe square roofs makes S{z) mulfi- 
valued. One is lead fo resolve fhis mulfi-valuedness, by considering m copies of fhe complex plane, 
one copy for each possible value of S{z) for given z. Lef us denote fhe m possible values of S{z) by 
5«(z), wifh i = 1,... ,m. Each copy of fhe complex plane musf also be provided wifh branch cufs, 
in order fo fully resolve fhe ambiguity. One should imagine cuffing fhe sheefs along fhe branch cufs. 
Afler cuffing fhe sheef, each branch cuf has fhe fopology of a closed curve. On copy i of fhe complex 
plane, fhe Schwarz funcfion is designated fhe value 5^*^ (z) for given z. Then fhe differenf copies of 
fhe complex plane are glued fogefher along fhe branch cufs, corresponding fo fhe facl fhaf as 2; crosses 
fhe branch cuf on sheef i fhe funcfion 5 ^*^ (2;) will coincide smoofhly wifh fhe value of function 5(2;) 
on some ofher sheef, j. 

For fhe example of fhe ellipse fhere are fwo sheds and 

^(1,2) . ) ^ (rf + rl)z ± 2rir2y/z2 + rj - rj 2) 

The square roof is faken in a way dicfafed by convention V re*^ = for —tt < 0 < tt. A 

branch cufs are drawn on bofh sheds belween (assuming ri > r 2 ). And fhe fwo sheds are 

glued along fhe branch cuf, fhe upper bank of fhe branch cuf on firsf shed fo fhe lower bank of fhe 
branch cuf on fhe lower shed. 

The algebraic Riemann surface. A, fhaf is produced by considering fhe complex curve F( 2 ;, 5'(2;)) = 
0, is also called an algebraic curve [16] . If is composed of m copies of fhe complex plane called 
sheds equipped wifh branch cufs. On each shed, away from fhe branch cufs, one defines a holo- 
morphic local paramder given by 2:. On fhe branch cufs, buf away from fhe branch poinfs, fhe local 
parameter spans more fhan one sheef. If zq is a simple branch poinf on shed i, fhen fhe local parameter 
around fhe poinf is ^/z — zq, spanning fhe fwo sheefs fhaf are connecfed af fhe branch poinf. 

We may define a funcfion on fhe algebraic curve A by specifying ifs values for any 2: and any shed 
i. If is fhus convenienf fo infroduce fhe nofafion z for a poinf on fhe A, where z = {z, i}, specifying 
fhe complex coordinafe, z, of fhe poinf, and fhe sheef, i. For example, S{z) is nafurally defined as: 

5(z) = 5({z,f}) = 5«(z). (3.3) 

Finally, note fhaf we shall called firsf sheef, i = 1, also fhe ’upper shed’, and all ofher sheefs, collec¬ 
tively as ’lower sheds’. 

The anli-analylic function r(z), defined as: 

r(z) = S{z), (3.4) 

is, in facf, a one-fo-one and onfo map from fhe Riemann surface fo ifself. This is due fo fhe facf fhaf 
if is ifs own inverse r(r(z)) = z, which is a consequence of unifarify (Eq. (2.14)). The sel of poinfs 
invarianf wifh respecf fo r is fhe inferface of fhe droplef, r(z) = z 5(z) = z. The anfi-holomorphic 
involufion, r, allows us fo equivalenfly describe fhe Riemann surface as a Schoflky double [17]. 

To describe fhe Schoflky double, consider fwo copies of fhaf parf of fhe complex plane fhaf lies 
oufside fhe droplef. Eel us denole Ihese fwo copies by and respedively. The Schoflky 
double, V is defined fo consisl of fhe fwo copies of fhe exlerior of fhe droplef, wifh local holo- 
morphic coordinates inheriled from fhe algebraic curve A, fhrough a one fo one mapping v, defined 
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Figure 2: The construction of the Riemann surface as an algebraic curve, A, is drawn on the left. The 
surface is three sheeted. The two banks of each of the branch cuts are numbered, such that banks that 
are glued together on different sheets are given the same number. A meandering closed path, through 
the Riemann surface is also drawn. The dotted line is the interface. The algebraic curve A is mapped 
onto the Schottky double V, in the middle drawing. The map of the meandering line is also drawn. 
The two connected components of the interface are now labelled A and B. On the right is a torus of 
genus 1 drawn as a doughnut. The interface with its two connected components, A and B, are drawn 
there, along with the meandering path. The torus is topologically equivalent both to the algebraic 
Riemann surface, A and the Schottky double V. 


below, from the algebraic curve A to the Schottky double V. The mapping v is defined as follows: 


( {z A} i = 1 and 2 outside C 
\ {r( 2 ;),|} otherwise 


(3.5) 


The upper sheet of the Schottky double, is usually called ‘the front side’, while the lower 
sheet, , is commonly known, apparently, as ‘the back side’. 

Note that since t{z) is anti-holomorphic, the holomorphic coordinate on is z, rather than z. 
Thus, by definition, an analytic on the Schottky double, is an analytic function of the variable z on 
and an analytic function of the variable z on The function must also be analytic across the 
interface. Another way to describe an analytic function, /, on the Schottky doubleis to say that / is 
analytic on the Schottky double, V, if 9 = f ° is analytic on the algebraic curve, A. 

As examples of analytic functions on the Schottky double, consider the function id{z) defined as 
fd({z,t}) = -2, id{{z,l.}) = z . Anofher example (wifh some abuse of nofafions) is fhe function 
S{z), defined as 5({2;,t}) = S{z) and 5({z,|}) = z, wifh apologies for using fhe same symbol 
S, for fhe Schoffky double funcfion, jusf defined, and fhe Schwarz funcfion. The fwo funcfions are 
disfinguished by fhe fhe facf fhaf fhe former fakes fhe combined nofafion {z, i} as an argumenf while 
fhe laffer lakes a complex number as fhe argumenf. In fhaf respecf, we shall offen use whaf we shall 
ferm as ‘fronl-side nofafions’. Namely, we shall denofe a function on the Schottky double by the 
same notation defined on fhe algebraic curve, if bofh funcfions agree on Thus, for example, 

fhe funcfion id above, would also be denofed by z. If fhe function is analytic fhere is no ambiguify 
infroduced in fhese nofafions, since fhere is a unique confinuafion of fhe function from lo 
on fhe Schoffky double jusf as fhere is a unique analyfic continuation from fhe upper sheel of fhe 
algebraic curve lo all lower sheels. 

Lei us make note of fhe non-frivial cycles on fhe differenl models of Riemann surfaces we have 
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Figure 3: The choice of cycles for a genus 2 Riemann surface. The b cycles draw a path, one part 
of which is on the front side and the other on the back side. The a cycles are drawn exactly on the 
contours connecting the front and back sides. 


defined. If the the interface has more than one connected component, the topology of the Riemann 
surface becomes non-trivial, in that there are cycles which cannot smoothly be contracted to a point. 
Two cycles are inequivalent if they cannot be deformed one onto the other smoothly. If the interface 
has g + I connected components, then we may define g such non-frivial inequivalenf cycles as g of 
fhe g + I connecfed componenfs. The lasf connecfed componenf is acfually equivalenf to the union 
of the first g, and thus may not be considered independent. These g cycles are denoted by aj, with 
j = 1,... ,g. In addition, g independent cycles, bj, may be defined, which infersecf fhe cycles aj, 
respecfively, wifh infersecfion number 1. This facf is wriffen formally as: 


di * bj 




(3.6) 


On fhe Schoffky double, fhe pafh of fhe cycle bj may be considered fo sfarf on fhe fronf side on fhe 
g + 1 connecfed componenf of fhe inferface, fhen confinues fo fhe back side infersecfing cycle aj 
finally connecfing fo fhe sfarting poinf fhrough fhe back side. The sifuafion is shown in Fig. 3. 

Finally lef use describe how fhe Hamilfonians, Hj., are defermined by fhe Riemann surface V 
(or equivalenfly A). The Hamilfonians, Hk, are mulfi-valued funcfions on fhe Riemann surface. Their 
derivafives wifh respecf fo z, are single-valued, due fo Eqs. (2.41), (2.29) and fhe single-valuedness 
of fhe Schwarz funcfion. Lef us define 



Hk T Hk zj{i) _ Hk Hk 

2 ’ ^ 


(3.7) 


which are fhe Hamilfonians conjugafe fo fhe real and imaginary parts of tk, respecfively. Denofe 
fhe real and imaginary parfs of tk by and , respecfively. We shall use fhe convention fhaf 
Ho = H/2, where by = 0, = H/2. 

The following are well-defined (single-valued) meromorphic differenfials on fhe Riemann surface: 


dH^ 


{R/I) 


dH, 


(R/I) 


dz 


-dz = 


dt 


dS 

{R/I) 


dz. 


(3.8) 


Due fo (2.44) and (2.41), fhe differnetial has a pole of order k + l&t infinify. Indeed: 


<'<’ = -^^^ + 0(1/."*), = (3.9) 
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where s = 2 ;“^ is the loeal parameter around infinity and the ± sign in the expansion of is 

to be taken respeetively with whether the expansion is around oO'|- or 00 ^. The differential of the 
Hamiltonians have no further singularities. This last faet ean be surmised from (3.8) and (2.19). 

A meromorphie differential with speeified singularities, sueh as preseribed by Eq. (3.9), is unique 
up to an addition of a holomorphie differential, of whieh there are g on a genus g Riemann surfaee. 
Sinee we ean add to any linear eombination of them with eomplex eoeffieients, the spaee of 

differentials obeying (3.9) is 2g dimensional over R. To fix eomplefely, one may explore fhe 

infegral over any eyele of sueh a differential. We firsf note fhaf Eq. (2.31) implies fhaf ^ S(z)dz over 
any eyele is purely imaginary (sef a = 6 in fhaf equafion). This fael, when eombined wifh (3.8) gives: 

= Re(f) =0. (3.10) 

fbj 

This eonsfifufes 2g real eonsfrainfs on whieh eomplefely fix fhe ambiguify in determining 

fhem for any k. This means fhaf dH^ can also be fixed uniquely given V or equivalenfly, A. As a 
eonsequenee we shall wrife in fhe following Hf^{A) in oeeasions where we shall wanf fo sfress fhaf 
fhe Hamiltonians are determined from the Riemann surfaee. 


Re(f) 

la, 


4 Dispersive Quantization and the two dimensional Toda lattice 

We shall now want to proeeed and quantize the Eaplaeian growth dynamies, eventually obtaining 
quantum field fheory, wifh fhe aim of using mefhods of infegrable and eonformal field fheory fo gain 
heifer underslanding of fhe fraelal paflerns in Eaplaeian growfh. The program of quanfizafion requires 
in fhis ease fwo sfeps. The firsf sfep is removing fhe dispersionless limif from fhe fwo dimensional 
redueed dispersionless Today hierarehy fhaf was presented above as fhe hierarehy defermining Eapla¬ 
eian growfh. This gives a elassieal nonlinear infegrable sysfem, fhe (dispersionful) fwo dimensional 
Toda laffiee. Eormally, fhe proeedure is very similar fo quanfizafion, Ihus if is dubbed here as ‘dis¬ 
persive quanfizafion’, allhough if eonneels fo elassieal problems. If is only Ihen fhaf we ean apply fhe 
seeond sfep, fhaf of physieal quanfizafion, fo oblain a quanlum field fheory. 


4.1 The two dimesional Toda Hierarchy 


The fwo dimensional Toda hierarehy [18], is an infegrable sysfem of nonlinear equations generalizing 
fhe more familiar Toda laffiee. The latter being a one dimensional sysfem of springs, having exponen- 
lial foree law F = Fq ■ — 1 )^ where F is fhe foree fhaf fhe siring exerls when if is exlended 

a dislanee Ax. The nonlinear equation is fhe Newlon equafion for fhe springs: 


m 


d^Xn 

dF 


Fo- 



(4.1) 


The fwo dimensional Toda laffiee, despile fhe name, does nol generalize fhe problem of exponenlial 
springs fo fwo dimensions, bul is ralher a formal generalizalion, in whieh time ilself is Irealed as a 
eomplex variable: 


dXn 

dtidii 


_ —1 ^n + 1 


(4.2) 


Here we have sef m, Fq and xq fo unify, as fhe physieal inlerprefalion in lerms of springs was already 
losl when we eomplexified lime, and so Ihere is no real use of keeping Iraek of physieal unils. 
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The equation of motion for Laplaeian growth, when treated as equations determining the evolution 
of a algebraie Riemann surfaee, are identieal to the equations that are obtained by eonsidering the slow 
time modulation of the two dimensional Toda lattiee hierarehy. Here one assumes that a solution to 
the two dimensional Toda lattiee has the property that it is rapidly oseillating in time (whieh itself a 
eomplex parameter, denoted here by ti), in sueh a way that the envelope of the oseillations varies on a 
mueh larger seale than the period of oseillations. If one parametrizes the envelope of oseillations using 
the eorreet variables, the equations for these parameters aequire an elegant form, whieh resemble a 
semielassieal limit of the full hierarehy. In the present ease these semielassieal equations are Eqs. 
(2.54), defining the Laplaeian growth evolution, while the equations for the two dimensional Toda 
hierarehy are obtained by replaeing the Hamiltonians Hk by differential operators and the Poisson 
braekets by commutators. This section, together with appendices A,B, C, is devoted to showing how 
this connection, described above in heuristic terms, may be derived rigorously. 

We must first define fhe whole fwo dimensional Toda lattice hierarchy, wifh ifs infinile number of 
complex limes, and fhe respeclive Hamillonians. The hierarchy may be defined [18] making use of a 
Lax operalor malrix, L. Lor fhe fwo dimensional Toda lattice, L is an infinile malrix depending on fhe 
limes, having fhe properly: 


Lnm = 0 for m > n + 1. 


(4.3) 


Each lime ti has associafed wifh if a Hamiltonian with matrix elements: 




{L )mn Tfl ^ n 
m = n 
0 m < n 


(4.4) 


here, e.g., {L^)mn is the mn element of the matrix L raised to the power 1. The time ti, has the 
as the associated Hamiltonian. Note the similarity of (4.4) with (2.34), (2.36). 

The nonlinear equations of the hierarchy may be obtained by demanding: 


dL^ 

dtk 




dL^ 

dik 




(4.5) 


As a consequence [19] of (4.5) and the definition of the zero-curvature conditions hold: 


dL(^) 

dtk dti 
dlik) 

dtk dii 


= 0 

= 0 . 


(4.6) 

(4.7) 


We shall be interested in a special reduction of the two dimensional Toda hierarchy. Namely, we shall 
assume the the Lax operator obeys the following condition: 

[L,Lt] = l. (4.8) 


In the dispersionless limit this equation becomes (2.21). 

To demonstrate how (4.6), (4.7), give rise to a whole hierarchy of nonlinear integrable partial 
differential equations, let us derive the first member of the hierarchy, Eq. (4.2), from them. Lirst note 
that any matrix satisfying (4.3) may be written as: 


1 

E “S’®®””. 

m=—oo 


(4.9) 
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where t aeting on a veetor whose elements are V’n gives a veetor whose elements are nV’n> while 
the operator when aeting on the same veetor gives a veetor whose elements are ilJn+m- 

= V’n, (iV)n = n'lpn, = V’m+n- (4.10) 


Note the similarity between Eq. (4.9) with (2.12). In faet the latter equation is the semielassieal 
limit of the former if we assume that in the limit p —)■ vp. The similarity between the dispersionless 
limit (the limit in whieh one deseribes the modulations of fast oseillating solutions) and a semielassieal 
limit, will present itself repeatedly below. The rigorous basis for the similarity will be treated in the 
next subseetions, supplemented by the appendiees. 

From Eqs. (4.9), (4.4) one may obtain as follows: 

+ (4.11) 


Setting A: = ( = 1 in (4.7) one obtains the equation: 


da, 


( 1 ) 


(f) , da^o \i) 


I I - 


m. 


dti 


+ 


dti 


dti 


paY\i)al^^\i) - aY’{i)a'^Q\i + 1) ai'^(f)a^'^(f) - a\"’{i)a''Q’{i + 1) . 

I t- “1“ C- . 




dti 


,(i) 


(i-l) 


.(1) 


(i) 


+ (4.12) 




For the operator on the left hand side to be equal to the operator on the right hand side, the free terms 
(the terms in the braekets on the left hand side and the first two terms on the right hand side) on the 
right and left hand sides must be equal, as well as the eoeffieient of and e“^, on both sides of the 
equation, respeetively. 

To obtain (4.2), one writes: 


,(i) 


(n) = e 


tOn 


^n+1 


(4.13) 


The the equation relating the eoeffieient of and e ^ on both sides of (4.12) lead to the same equa¬ 
tion: 

-z 2 ^ + -^{Xn - Xn+i) = ao^^(ra) - ao\n + 1 ) 

Oti uti 

dao\n) da‘il\n) ^ 
dii ^ dti 

Adding the ti derivative of (4.14) to its eomplex eonjugate, one obtains: 
d^ 

r, r,- (Xn - Xn+l) = ^n_^^x„+i x„+2_2Q^n Xu+1^ 

utiOti 

whieh is just what is obtained by subtraeting Eq. (4.2) from the same equation when n is shifted up 
by 1, (that is n —)■ n + 1). Namely, if we assume that —)■ 0 fast enough as n —)■ ±oo, then Eq. 
(4.16) is equivalent to (4.2), sinee in that ease we ean use, e.g., Xn = J2m=n “ Xm+i- 


(4.14) 

(4.15) 
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4.2 Multi-Periodic Solutions 

As mentioned above, the redueed dispersionless two dimensional Toda lattiee hierarehy, whieh de¬ 
termines the evolution Laplaeian growth ean be obtained by taking the dispersionful redueed two 
dimensional Toda lattiee hierarehy and eonsidering slowly modulated waves. In order to be able to 
diseuss these modulated waves, we must introduee the unmodulated periodie solutions. In partieular 
we must diseuss how these are eonstrueted. To do this it is useful to make use of the Baker Akhiezer 
funetion [20]. 

The Baker Akhiezer funetion, is an eigenfunetion of the Lax operator, with eomplex eigenvalue z: 

L'lp = z'lp. (4.17) 


The time dependenee of V’ is defined naturally using the Hamiltonians: 

^ A: = (4.18) 

Otk otk 

For eaeh z there may be several eigenfunetions of L. Suppose that, generieally, there are m of 
these eigenfunetions for eaeh z. We may take m eopies of the eomplex plane C and assign to eaeh 

eopy one of the eigenfunetions. It turns out that there is a large elass of L’s sueh that by drawing 

braneh euts on the different sheets and gluing the sheets together along the braneh euts, the funetion V’ 
becomes a smooth function of z, except at infinity, and at a set of discrete poles. In fact, Krichever [21] 
has shown that an effective way to construct multi-periodic solutions is to start with an algebraic 
Riemann surface, and define fhe Baker-Akhiezer funcfion according fo some analyfic condifions (see 
also Ref. [22] for an overview of fhis subjecf). The Lax mafrix L, can fhen be consfrucfed from fhe 
knowledge of fhe Baker-Akhiezer funcfion. Since fhe elemenfs of fhe mafrix L are fhe nonlinear 
fields of of fhe fwo dimensional Toda hierarchy, fhe lax mafrix also encodes solufions fo fhe nonlinear 
infegrable equafion. 

If fhe algebraic Riemann surface in secfion admifs an anfi-holomorphic, r, fhen fhe Riemann 
surface can be given fhe sfrucfure of a Schoflky double. Anficipafing fhe connection fo fhe Riemann 
surfaces we encounfered in Laplaeian growfh, we assume fhe existence of such a anfi-holomorphic 
involufion and fhe problem of consfrucfing mulfi-periodic solufions fhus becomes fhe problem of 
consfrucfing Baker-Akhiezer funclions on a given Schoflky double. 

The defails of how fo define a Baker-Akhiezer funcfion from a sel of analytic properfies is given 
in appendix A. The facl fhaf such a funcfion fhen solves (4.18) is fhen shown in appendix B. The 
informalion perfinenf fo fhe sequel is fhaf for each Schoflky double a procedure yields a mulfi-periodic 
solution of fhe form: 


Ur 


i{to) = /) 





(4.19) 


where u)k are g dimensional complex veclors and fm are funclions from fo C, which are periodic 
wilh periods e* and Bei'. 


fm{v -h ei) = fm{v), fm{v + BCi) = fm{v). (4.20) 

Here e* is fhe unil veclor in fhe f-lh direclion, {ei)j = 6ij, and B isagby g mafrix. Many more defails 
are given in appendices A and B. However, fhe main message of Eq.(4.19) is fhaf a complicated mulfi- 
periodic solulion of fhis form may be oblained from a given Schoflky double. In facl, fm is given as 
a combinafion of Riemann Ihefa funclions associaled wilh fhe Riemann surface. 
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4.3 Modulations Equations 

Given a multi-periodic solution, corresponding to some Schottky double, or equivalently, to a al¬ 
gebraic Riemann surface, one may consider the effect of applying any perturbation to the nonlinear 
equations. The effect of such a perturbation would be to slowly change the nature of the multi-periodic 
solution. Namely, the amplitude, phase, average value and frequency would change. After a while, 
the solution would be described by a different Schottky double, than the one describing the solution at 
initial times. We may describe the modulation of the multi-periodic wave as dynamics of the Schottky 
double or the algebraic curve. 

The form of the dynamics of the Schottky double naturally depends on perturbation applies. It 
turns out, however, that if the strength of the perturbation is sent to zero, the original, unperturbed, 
solution is not recovered. In other words, the limit of vanishing perturbation is singular. Nevertheless, 
the limit of zero perturbation is universal, in that for sufficiently well behaved perturbations, taking 
the amplitude of the perturbation to zero yields the same universal Whitham dynamics on Schottky 
doubles. This leads to the interesting situation in which there is a natural dynamics defined on Schof- 
fky doubles or, equivalenfly, algebraic curves. For fhe case of fhe reduced fwo dimensional Toda 
lattice, fhe universal dynamics obfained is fhaf of Laplacian growfh, which was already described as 
fhe dynamics of fhe algebraic Riemann surface A. 

Whifham [23-27] developed some of fhe firsl fools fo describe such universal modulafed dynam¬ 
ics for nonlinear wave. Flaschka, Fores! and McLaughlin connecfed fhis evolufion fo concepfs in 
algebraic geomefry, while Krichever [28] was able fo derive fhe dynamics rigorously making use of 
an averaging mefhod fhaf he had devised. The resulfing dynamics has an infegrable sfrucfure (as we 
have seen in Laplacian growfh above). Such dynamics are somefimes called dispersionless hierarchies 
or Whifham hierarchies. Appendix C is devofed fo a review of Krichever’s work [28], following fhe 
more accessible presenfafion of Ref. [23,27]. 

The modulafion equations obfained by fhis procedure are Eqs. (2.51). They are obfained by es¬ 
sentially averaging fhe zero curvafure condifions, Eqs. (4.6), (4.7), wifh fhe Baker-Akhiezer funclion, 
and ifs adjoin! over many periods of fhe mulfi-periodic solufion, buf on a scale much smaller fhan fhe 
fypical scale of modulafion. The defails are given in appendix C. 

Eel us nofe fhaf Eqs. (2.51) can be solved by faking a Riemann surface which is limes independenl. 
Indeed, as was describe in seclion 3, fhe Hamiltonians Hk can be found by considering fhaf dH^ is 
a unique differential on fhe Riemann surface, such fhaf if fhe Riemann surface does nol change Ihen 
bolh lefl and righl hand side of (2.51) are zero. Neverlheless, such a solution is obviously nol a 
solution of fhe Eaplacian growfh problem. Malhemalically, we see fhis by nofing fhaf we are seeking 
solutions in which Ihere is a Schwarz funclion fhaf generates fhis Hamilfonian (Eq. (2.41)) and obeys 
fhe unilarily condition (Eq. (2.14)). The facl fhaf given initial conditions fwo solulions exisl, is relaled 
fo fhe facl fhaf fhe Whifham equations are a singular limil obfained by adding a perlurbalion fo fhe 
equations and sending fhis perlurbalion fo zero. The solufion which corresponds fo Eaplacian growfh 
is fhe non-lrivial solution in which fhe Riemann surface is lime dependenl and in which fhe unilarily 
condilion is preserved by fhe evolution. 

5 The Classical Korteweg-de Vries Limit 

Afler discussing how fo remove fhe dispersionless limil from fhe infegrable sfrucfure of Eaplacian 
growfh to oblain fhe (dispersionful) reduced fwo dimensional Toda lattice hierarchy, we would like fo 
quanlize fhe latter. A lol more is known, Ihough, on fhe quanlizalion of Korleweg-de Vries equation 
fhan fhe quanlizalion of fhe fwo dimensional Toda lattice. In facl, quanlum Korfeweg-de Vries is fhe 
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integrable structure underlying conformal field theory, which makes this quantum integrable system 
very appealing, due to the rich structure and the applicability of results of the representation theory of 
the Virasoro algebra. 

Due to these consideration, we opt to first treat a limit of Laplacian growth, which has the inte¬ 
grable structure of Korteweg-de Vries associated with it [10]. Only then will we consider quantization. 
The advantage of this approach is that we will be able to use results from conformal field theory after 
quantization, the drawback is that it is not clear to what extent the fractal nature of the growth extends 
to the Korteweg-de Vries limit. 

5.1 Narrow Finger Limit of Hele-Shaw Flows 

If one focuses on a region around a Hele-Shaw finger, the integrable structure of the two dimensional 
Toda lattice simplifies in the limit of a very thin and long finger. Let us orient the x axis to point in 
the long direction of the finger (see Fig. 4). The Korteweg-de Vries limit occurs when the finger is 
symmetric to reflection across the x-axis and when the curve is described at intermediately large and 
negative x as 


J-3/2 



(5.1) 


-h .. . , 


where xq 3> yo, so that the finger is long and narrow. Eq. (5.1) describes the shape of the finger 
at X ~ xo, and may have more complicated features at x <C xq. Nevertheless, it is assumed that 
throughout the region, for all points on the interface, y <C x. Due to this /^(^<y7) is approximately real 
in the entire region. 



Figure 4: A small part of the droplet is magnified, to reveal a long and narrow, possibly multi- 
connected, finger, symmetric with respect to reflections about the real axis. 

The Korteweg-de Vries description of the finger holds for any rapid evolution of the interface at 
|x| <C Xo, while the asymptotic behavior at x ~ xq, Eq. (5.1), does not change. Namely, in Eq. 
(5.1), we may assume that xq and yo are time independent. Furthermore, since the Korteweg-de Vries 
description only holds for the confined region x <C xq, we may assume that Eq. (5.1) describes the 
asymptotic behavior of the not only at x ~ xq but all the way to x —)• — oo. 

To describe such a finger, we first define a re-scaled and shifted and t given by (p and i, respec¬ 
tively. Namely, we define (p{<p) = aip + b, i = at + j5, where the way a, b, a, /3 are chosen will be 
described below. Eet: 



(5.2) 
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We may now expand hi as a Laurent series in (p around infinity, sueh that the asymptoties deseribed 
in Eq. (5.1) is reeovered.The simplest ehoiee is to take: 

k = ‘f>^+ ( 5 - 3 ) 

j=-oo 

and 

h-t = O • (5.4) 

The ehoiee of a and b in (p{ip) = aip + b, is made sueh that the free term in (5.3) is missing and the 
highest order term has unit pre-faetor. The faet that li has only even powers in an expansion in ip and 
only odd powers, eorresponds to the symmetry of the finger wifh respeef fo refleefions aeross fhe real 
axis. Indeed, fogefher wifh (5.4), fhe refleelion symmefry is effeeled by (^ —)■ —(p, whereby f^ —)■ Ip 
while hi —)■ —hp The fael fhaf Ibis way fo respeef fhe refieefion symmefry is aefually implemenfed by 
fhe eonformal map ft in fhe asympfofie region, is nof guaranfeed. Neverfheless, if may be expeefed fo 
hold fairly generally, when fhe shape of fhe inlerfaee away from fhe finger does nof induee a sfrongly 
asymmefrie pressure field on fhe finger. We shall nof dwell on Ibis poinf, as we are mainly inferesfed 
in fhe exisfenee of fhe Korfeweg-de Vries limif under some eondifions, rafher fhan demanding fhaf if 
hold under generie eondifions. 

Insfead of fhe Sehwarz funefion if is more appropriafe fo define fhe following eombinafion: 

y{x) = hi{ipi{x)), (5.5) 


where ipi{x) is fhe inverse funefion of If 


hiViix)) = X. (5.6) 

The funefion y{x) is fhe analyfie eonfinuafion in x—plane of a funefion whieh gives fhe y eoor- 
dinafe of fhe finger as a funefion of fhe x eoordinafe, for x fhaf belongs fo fhe droplef. Sueh x’s fhaf 

below fo fhe finger form segmenfs on fhe real axis. These segmenfs are braneh eufs for fhe funefion 

y{x). If one approaehes fhe braneh euf from above or below one obfains fhe ±y(x), respeefively. 
Namely as x approaehes fhe real axis we have y{x) = —y{x). Sinee y is real we ean lake fhe eomplex 
eonjugafe of fhe leff hand side of Ibis equalion and oblain: 

y{x) = -y{x). (5.7) 

This equalion holding for x approaehing fhe real axis. 

Around fhe eonlour we have osfensibly y{x) = Xln o^nWx — Aj)”. To obey (5.7), a 2 n have fo 
be purely imaginary, bul y is real on fhe real axis and so a 2 n = 0. Thus, y{x) has only odd powers 
of yjx — \i whieh implies fhaf y‘^{x) is regular, as a funefion of x, everywhere in fhe eomplex plane. 
Namely, if is an entire funefion. In addilion, y‘^{x) ~ x‘^'^~^ as x —)• oo, whieh, fogefher wifh fhe fael 
fhaf if is entire, suggesls fhaf if is a polynomial. Considering fhaf fhe zeros of Ibis polynomial are fhe 
braneh poinls, one obfains: 


25+1 

y^{x) = n 

i=l 


(5.8) 
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This equation describes a Riemann surface, through a complex curve. Curves which are defined 
through = Pn{x), where is a polynomial, are called hyper-elliptic. Thus, in the Korteweg-de 
Vries limit, the Riemann surface is a hyper-elliptic curve. 

Given the definition, Eq. (5.2), of and h^, Eq. (2.21) implies the following Poisson bracket for 
these objects: 


{h,hi} = h (5.9) 

where here the choice of a in f = at -I- /3 is made as to have 1 on the right hand side. The choice of /3 
remains arbitrary. Eqs. (5.3), (5.4) leads to the following expansion for y: 

J 

y{x) = Y, + 0(z-3/2). (5.10) 

k=l 

with ii = i. All the f^ 2 fc+i are real parameters, as the symmetry of the finger dictates. 

In a manner following closely to the way in which Eq. (2.29) was derived from (2.21), one may 
derive from (5.9) the following equation: 


Define now the Hamiltonians 


dy{x) dtfi{x) 
di dx 


H2k+1 — 




(5.11) 


(5.12) 


where the definition of , for any function g{(p) = J2T=-oo is as follows: 

OO 

= f+ E“‘A 

k=l 


Assume 


for k odd. Then, 


dik 


{xi,Hk}, 


dhj 

dik 


{ , Hk }) 


dyjx) 

dik 


dli{(p) d(pt{x) 

dik dik 


^=(pt(x) 


= {IhHk] 


( dit{(p) 

\ 


^ dlt{‘f>) diti^) 
d(f dik 




dHk , dit{(p) 
dip dip 


(p=pt{x) 


dHk 

dx 


(5.13) 

(5.14) 


(5.15) 


(5.16) 


(5.17) 


This equation is an evolution equation in time tk for y, and has the property that it retains the form 
(5.10). This shows that the times ik may be viewed as the analogue of the harmonic moments in 
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the narrow finger limit. First, eaeh tk is eonserved when any on of the other times is ehanged and 
seeondly, the ik are quantities eharaeterizing the shape of the interfaee. Indeed, the funetion y{x), for 
whieh the ik are Laurent eoeffieients, ean be found by direetly from the shape of the interfaee, as it is 
the analytie eontinuation in the eomplex x—plane of the funetion whieh gives the y eoordinate of the 
interfaee given the x eoordinate. 

5.2 Dispersive Quantization of Korteweg-de Vries 

Just as we have applied ‘dispersive quantization’ the the dispersionless limit of the two dimensional 
Toda lattiee equations, we shall now be interested in applying the same proeedure the the dispersion¬ 
less Korteweg-de Vries equations (5.15). We shall at this point do away with the tildes over the time, 
over the (p variable, and the Hamiltonians, sinee we shall not treat the two dimensionless Toda any¬ 
more, and as sueh there should be no eonfusion about whether the objeets diseussed belong to the 
Korteweg-de Vries or the two dimensional Toda hierarehies. 

The Korteweg-de Vries hierarehy may be defined (See Ref. [29] and referenees therein) by eon- 
sidering the Lax operator, L, a differential operator in ti, given by: 

^ = + (5-18) 

One may write the square root of L as a formal power series in dt^: 

OO 

+ bm{ti)d^/^, (5.19) 

m=0 

where 0^“”^ is the formal inverse of dt^ raised to the power m. The funetion bm{ti) can be found, 
order by order, by demanding = L. 

A set of Hamiltonians are then given by 

= (5.20) 

where the (L^+^/^)_i_ denotes only positive powers of dt^ in the expansion of whieh itself 

eontains both positive and negative powers. The equations for the hierarehy follow from the Lax 
equations: 

dt,L=[UHk]l (5.21) 

The dispersionless limit of this problem beeomes the problem that has been diseussed in the previ¬ 
ous subseetion. To show this one may follow mueh the same as deseribed in for the two dimensional 
Toda lattiee. In faet the result may be obtained by taking an appropriate limit of that proeedure. 
Namely, the long and narrow linger limit, in whieh the Riemann surfaee beeomes hyperelliptie, and 
only real odd times deseribe the evolution. As sueh, we shall not repeat those steps. We note, that the 
original literature eited here on the matter [28,30,31] either deals direetly with the Korteweg-de Vries 
problem, or with elose relatives, whieh are easily redueed to the Korteweg-de Vries problem - mainly 
the Kadomtsev-Petviashvili hierarehy. 

The first member of the hierarehy gives the tautology dt^u = Ot^u.The first non-trivial member 
of the hierarehy is given by: 

^dt^u = + udt^u. (5.22) 
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Alternative to the Lax eonstruetion, this equation may be obtained by defining the following Poisson 
braekets: 


{uix),u{y)} = 2(u(x) + u(y))6'(x - y) + S"'(x - y), (5.23) 

and allowing for evolution with the Hamiltonian: 

= u^{ti)dti, (5.24) 

upon whieh, Eq. (5.22) beeomes equivalent to: 

dt,u = {u,H 3 }. (5.25) 

The tautologieal equation dt^u = dt^u is generated by Hi given by 

Hi = J u{ti)dti. (5.26) 

All higher order flows defined by (5.21) may be obtained using the Poisson braeket (5.23) by defining 
appropriate Hamiltonians (See Ref. [29] and referenees therein). 


5.3 Inverse Scattering 

The method of inverse scattering was first invented for the Korteweg-de Vries equation [32] and later 
expanded to become a wide-ranging subject (see, e.g., Ref. [29] and references therein). It is useful 
to review briefly the inverse scattering method, since it provides a link between the classical multi- 
periodic solutions to the quantum states one obtains in solving the integrable quantum Korteweg-de 
Vries equation. 

To present the method of inverse scattering, first let us consider the periodic problem. Namely, we 
consider the problem in which the field u{ti) is periodic with period 27r. Next, we consider the Miura 
transformation: 


\dti) dtj' 


(5.27) 


The Poisson bracket for the field (p is then given by: 

{(j){ti),(j){t'i)} = ^sign(ti -t'l), 


(5.28) 


and can be shown to be equivalent to (5.23). Namely, substituting (5.27) into (5.28) one obtains (5.23) 
up to full derivatives. The latter are determined by requiring that the Poisson bracket of <p must be 
obey the Jacobi identity. 

The Korteweg-de Vries equation, Eq. (5.22), translates into the following equation for p: 


dp _ d^p ( Of/* \ ^ 

dtz dt\ V9fiy 


(5.29) 


which is a version of the modified Korteweg-de Vries equation (the standard modified Korteweg-de 
Vries equation is written for the field v = pi). 
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The equation for (j), Eq. (5.29) ean be obtained for a zero eurvature eondition for a Lax pair. 
Consider the operators 


Ax = + \/A(T2;, (5.30) 

Bx = \/A(4A - 20'2)f7^ + + (j)'{4X - 2^'^))a, - 2^0"\/AcJy, (5.31) 


where the prime denotes derivative with respeet to ti. Then the equation 


A 

dti 




= 0 , 


(5.32) 


becomes equivalent to (5.29) as can be ascertained by computing the commutator explicitly. 

Eq. (5.32), which is equivalent to the dynamic equation of cj), Eq. (5.29), has exactly the required 
form to serve as the consistency condition required for there to be a solution, to the equations: 




(5.33) 

(5.34) 


where V’ is a two dimensional vector. The significance of this fact is that one may approach the 
problem of finding solufions of fhe modified Korfeweg-de Vries equafion, Eq. (5.29), or equivalenfly 
fhe Korfeweg-de Vries equafion, Eq. (5.22), as fhe problem of consfrucfing solufions fo (5.33) and 

(5.34) as Baker Akhiezer funcfions as done above. In fhis secfion we shall show how fhe Baker- 
Akhiezer function is related fo fhe mefhod of separafion of variables. 

The solution of (5.33) is in facl simple. We may write: 


At) 




/o' 


m, 


(5.35) 


where ^ denotes a pafh-ordered producf. One is lead fo define fhe shiff operafor: 


Mx{t'i,ti) = e-^*i 


(5.36) 


The problem of finding Bloch wave-functions V’(fi + 27r) = satisfying (5.33) and 

(5.34) is fhen fhe problem of finding eigenvectors of fhe mafrix Mx{2tt, 0). The objecf is called 
fhe Bloch multiplier of fhe Bloch wave-function V'(fi). Note fhaf fhe firsl elemenf of y), denoted by 
V'l, is also Bloch wave-functions for fhe following equation: 

+ tt(fi))V’i = At = AV'i, (5.37) 


an equation which is easily derived making use of Eqs. (5.33) and (5.27). A Bloch eigenfunction of fhe 
Lax operator is fhe Baker-Akhiezer funclion. Indeed, fhe Baker-Akhiezer funcfion is a solution fo fhe 
specfral problem, where fhe coefficienls appearing in fhe Lax operafor are mulfi-periodic funcfions. 

A special case in which fhe eigenvectors of Mx{27r, 0) are apparenf presenfs ilself for values of A 
af which Mx{27r, 0) happens fo be diagonal. The fhe Bloch wave-functions are given by fhe sfandard 
basis 


1 

0 



(5.38) 
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up to a times dependent faetor. We see that the first element of one of the Bloeh funetion vanishes. 
Namely the Baker-Akhiezer funetion vanishes at this A on one of the sheets. The matrix Mx depends 
on the times and as sueh the points A at which this matrix is diagonal also depend on time. Suppose 
that there are g such points. Indeed, on a Riemann surface of genus g the Baker-Akhiezer function has 
g zeros (see appendix A). These points form a set of g dynamical variables denoted by 7 *. Another set 
of dynamical parameter may be taken as the Bloch multiplier associated with the non-vanishing Baker- 
Akhiezer function at 7 *. Namely we take the upper left element, , as a set of g additional 

dynamical variables. As shown in Ref. [33] and reviewed below, these parameters have a Poisson 
structure that allows one to treat them as separated variables. We review the computation of these 
Poisson brackets in the following, by computing the Poisson brackets of the elements of M\{2'k, 0) 
for any A and then specializing to the points where Mx{2tt, 0) is diagonal. 

We wish to find fhe Poisson brackefs of fhe differenf elemenfs of Mx{27r, 0). Since we are infer- 
esfed in finding fhe Poisson brackef of any fwo mafrix elemenfs of Mx{2tt, 0), if is useful fo consider 
fhe simulfaneous operation of faking fhe fensor producf and fhe Poisson brackef. Namely for any fwo 
mafrices, a and b, one may consider fhe objecf {a®b}, which is an objecf which lives in fhe same 
space as a (g) b, and has mafrix elemenfs given by: 

(5.39) 

jusf as (fl (g) 

From now on we shall denote 

Mr = Mr i2Tr,0) (5.40) 

Wifh Ibis definilion we compufe: 

{Mx^Mfj,} = jj dxdx'Mx{2Tr,x) 0 Mfj_{2Tr,x'){Ax{x)^A^{x')}Mx{x,0) 0 Mfj,{x',0) 

(5.41) 

This equafion is a consequence of fhe definilion of M in terms of A (Eq. (5.36)) only, and may be 
derived by writing M in terms of a path ordered product: 

1 + / Ax], (5.42) 

%(-i) J 

and computing the Poisson brackets that appear in (5.41). 

An explicit calculation gives 

{AA(x)f’A^(x')} = 6'{x - x')az (g) (5.43) 

Inserting this into (5.41), writing 5'{x — y) = — y) and integrating by parts, one obtains: 

{Ma®M^} = (5.44) 

= J dxMx{27r,x) 0 M^{2 tt,x) az 0 (Tz, ^ ^ ^ l(g)A^(x) ® M^j,{x,Q), 

where here we used the fact that Mr, obeys similar equations to those that y) obeys: 

dxMr{x,x) = Ar{x)Mr{x,x'), dx'Mr{x, x) = —Mr{x, x') Ar{x), (5.45) 
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whereupon one may substitute /x or A for r. 

Finally, we need the following identity, which may be obtained by an explicit calculation: 


kz (8) (Tz, Ax{x) (g) 1 - l(8)A^(x)] 


(8) 1 + l(8)A^(x) , 


where 

„ ''■ + T~^ (Xz® 2 , 

-Rr = -f-X-1- 7 ( 0 -+ (8) CT_ + CJ_ (g) a+). 

T — T~^ 2 r — T~^ 

Substituting (5.46) into (5.44) and making use Eq. (5.45), one may write: 


= dxdx{Mx{2TT,x) ® ® 

R mJ , 


(5.46) 


(5.47) 


(5.48) 


This is the desired result. In the present, classical, terms, the equation encodes the Poisson brackets 
between the elements of M at different spectral parameters in terms of bilinears of M itself. The 
quantization of these relations will allow to compute the commutations relations for such elements in 
the same manner. 

An almost immediate consequence of (5.48) is that the trace of Mx, denoted as T^: 


Tx = tiMx, 


(5.49) 


Poisson commutes 


{Ta,T^} = 0, (5.50) 

an equation which is immediately obtained from (5.48) by noting that tr [{afb}] = {tr[a],tr[ 6 ]}, 
just as tr [a (g) 6 ] = tr[a]tr[ 6 ]. 

The relation (5.50) encodes the infinite number of conserved quantities in involution. Indeed, it 
turns out that the expansion of Tx around A —)■ 00 , yields the conserved quantities as the coefficients 
in the Laurent expansion. This proves that they are in involution. 

It is customary to define fhe following mafrix elemenfs of fhe mafrix Mr'. 


Mr = 


ttr 

Cj- 



(5.51) 


To obfain separated variables one firsf nofes fhaf from (5.48) fhe following relafion may be obfained 
by faking fhe appropriafe mafrix elemenf: 


{bx,b^} = 0. (5.52) 

Consider fhe poinfs af which Mr is diagonal. Af fhese poinfs, denoted by 7 * we have b.y- = 0. The 
following is fhen a consequence of (5.52): 


= 0 . 


(5.53) 
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We have found a set of variables whieh are in involution. As further dynamieal variables one may 
eonsider A, = a^.. Referring back to (5.48) and taking the appropriate matrix elements one concludes: 

{Ai,A,} = 0, (5.54) 

and 

(5.55) 

where b'^. = ^ I,—.y • Eq. (5.55) may be obtained if we postulate the following Poisson brackets: 

{Aj, 7 j} = ASijAiji. (5.56) 

We thus have: 

{A*, Aj} = { 7 i, 7 j} = 0, {logAj,log 7 j} = (5.57) 

The variables log( 7 i) and log(Aj) are thus canonical. To show that these are also separated vari¬ 
ables, one needs to show that each pair ( 7 ^, A*) traces out a one dimensional loop, whose shape is 
independent of all the rest of the canonical variables. Namely, one must show that Aj is a function of 
7 j, the function depending only on the conserved quantities. The existence of such a representation is 
a consequence of the fact that Aj is the Bloch multiplier associated with the Baker-Akhiezer function 
at 7 j. An explicit expression for the Baker-Akhiezer function is given in appendix A, from which 
A( 2 :) may be found. We thus have A* = A( 7 j), as required. 

6 Quantum Korteweg-de Vries 

Having reviewed the classical inverse scattering method, we are able to present the quantum version of 
which, discovered in Refrs. [11-13] based on results of Refrs. [14,15]. More context on the quantum 
inverse scattering can be found in Ref. [34] and is also treated briefly in Ref. [29]. We then go on to 
review the quantization of the Korteweg-de Vries problem based on quantizing the field u{ti), which 
leads to the connection to conformal field theory. An object which bridges the two approaches is given 
by the form factors, which has a central role in the proposed connection to Laplacian growth. 

6.1 Quantum Separation of Variables 

The classical separation of variables method for the Korteweg-de Vries problem was achieved by 
considering the monodromy matrix. Quantum Korteweg-de Vries can be defined by quantizing in the 
separate variables, in a method suggested by Sklyanin [35]. On defines the quantum analogues of 7 * 
and Aj as a set of operator 7 *, Aj, with the following commutation relations: 

[Aj, Aj] = [ 7 j, 7 j] = 0, Aijj = (6.1) 

To see that indeed (6.1) is the classical limit of the last equation of (5.57), we write the last equation 
in ( 6 . 1 ) as: 

glog(Ai)+log( 7 fc)-l-|[log(Ai),log( 7 fc)]-|-... _ gi7r;j(5ij-|-log(Ai)-|-log(%)-l[log(Ai),log(7fc)]-|-...^ |-g_ 2 ) 
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obtained by using Baker-Campbell-Hausdorff Lemma. Taking the logarithm of both sides, one ob¬ 
tains: 


[log Aj, log 7 j] = SijiTTh + — (6.3) 

The semiclassical limit of this equation is indeed (5.57), as desired. 

One may find a representation of this algebra by postulating some Hilbert space spanned by the 
joint eigenvectors of all the y^’s, which we shall denote by I 7 ), with 7 i| 7 ) = 7 i| 7 ). The operator Aj 
can then be represented on such states as: 

Aj I 7 ) = I 7 ). (6.4) 

To make contact with a more usual field fheory descripfion of quanfum Korfeweg-de Vries, namely, 
a description in which fhe field u{x) is quanfized, fhe variables Aj, 7 *, musf be written fhrough fhe 
quanfum field u{x) . In facl if is more advanfageous fo sfarf wifh fhe field Miura field, (j) (see Eq. 
(5.27)), and quantize if, and fhen perform a quanfum Miura fransformafion fo fhe field u. To fhis 
aim, we define an operafor analyfic in fhe upper half x-plane, which satisfies fhe commufafion 
relations: 

= ( 6 . 5 , 

The latter being an analogue of (5.28). 

In analogy wifh (5.36), one may define an operafor valued mafrix wave-function, ^A(a^) depending 
on fhe specfral paramefer A, as follows: 

^^(x) = exp 

where d±, dz are mafrices salisfying: 

giFuJz _ g— 

[cr„cr±] = ±2d±, [d+,cr_] = , (6.7) 

which in fhe limif /i —0 yields fhe usual Pauli mafrices. ^a( 3:) is a operafor-valued mafrix solution 
fo a quanfum analogue of Eqs. (5.33), (5.34). 

The rafher unusual way in which fhe paramefer h enfers info ( 6 . 6 ) fhrough (6.7) is necessary in 
order for 'Pa (a;) to produce, through the monodromy matrix, a set of quantum variables satisfying 
(6.1). Indeed, define fhe quanfum monodromy mafrix, M\ as fhe mafrix 'kA(®) af x = 27r: 

Ma = 4'A(27r). (6.8) 


e-2'^(^')cr+ + e2'^(^')d_) dx' 


(6.6) 


If can be shown [13], using fhe commufafion relations (6.5), fhe fhe following relation holds for M: 


wifh 


R (a/Vm) 0 1)(1 ® Mfz) = (1 ® ® 1)7? (a/Vm) • 


m 


( 

V 


lA - qX- 
0 
0 
0 


1 


0 0 ® \ 

A — A“^ q~^ — q 0 

q~^ — q A — A“^ 0 

0 0 q-^X-qX-^ / 


(6.9) 


( 6 . 10 ) 
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This relation is the quantum analogue of (5.48). Indeed, we set q = e*^ and write: 

= 1 - zTi (r(A) + + 0{h^) (6-11) 

If one substitutes (6.11) into (6.9) and uses the usual semielassieal expansion of the produet of oper¬ 
ators, AB = AB + '^{A,B} + ... one obtains (5.48) to first order in the expansion. Namely, the 
eommutation relations between the matrix elements of M, whieh are eneoded in (6.9) ean be viewed 
as a quantization of the eorresponding Poisson braekets, whieh are eneoded in (5.48). Thus, one may 
proeeed with a quantum version of the separation of variables method, based on (6.9). 

The proof of (6.9) is rather eumbersome and Ref. [13] (the third from the series of papers [11-13]) 
foeuses on this aspeet of the quantum Korteweg-de Vries integrability, the ealeulations are based on 
insights gained from Ref. [14,15]. 

Let us write: 


( A{\) B(X) \ 

- ( C{\) D(A) ) ■ 

We are now ready to write 7 * and A* through (j), albeit rather indireetly, sinee we shall write these in 
terms of A, B, C and V, whieh ean, in turn, be written through (j), due ( 6 . 6 ), ( 6 . 8 ). The following 
expansion of the operators may be shown: 

/ ,2\ °° °° 

B{X)=B{0)Ylil-^U A{X) = Y,X^^An, V{X) = Y,X^^Vn, (6.13) 
j V J n=0 n=0 

whieh defines ji, while A*, is defined by A{^i) = Aj. In fhe subslifufions A{'^i), fhe operator if'" is 
fo be plaeed fo fhe lefl of An respeefively. The eommufafion relafions (6.1) follow from fhe definition 
of 7 i, Aj, and from Eq. (6.9). 


6.2 Conserved Quantities 

We now pass to a deseripfion of quanfum Korfeweg-de Vries fhaf sfresses fhe quanfizafion of fhe field 
u{x). As fhe quanfum analogue of fhe Miura Iransformafion (5.27) one may lake: 

f{x) = 4'^{x)-. + (1 - 0 0"(x) - (6.14) 

The field ^T{x) is fhe quanfum analogue of fhe Korfeweg-de Vries field u. If is designafed as T sinee, 
as defined, if has fhe eommufafion relafions of fhe slress energy tensor of Iwo dimensional eonformal 
field fheory. Indeed, fhe eommufafion relafions for (p, Eq. (6.5), beeome fhe following eommufafion 
relafions for fhe field T: 

i[f{x),f{y)] = - (t{x) + f{y)) 2tt8'{x - y) + ^5'''{x - y), (6.15) 

whieh are bofh familiar from eonformal field fheory and eonslilule a quanfizafion of fhe eommufafion 
relations (5.23), where 

c=U-q{^ + ^ . (6.16) 
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The commutation relations (6.15) becomes the commutation relation (5.23) in the limit /i —0 if we 
assume that in this limit 

f{x)^^u{x), [•,•]-> •}. (6.17) 

Given Eq. (6.15), it is possible to construct an infinite set of mutually commuting quantities on 
the cylinder, 



(6.18) 

(6.19) 

( 6 . 20 ) 


where the normal ordering is defines as: 

: AB : (w) = ^ (f -A^A{z)B{w), 

2m J z — w 


( 6 . 21 ) 


and fhe infegral is over a small circle surrounding w. 

Thai fhese are mufually commufing quantifies may be ascerfained [36] by finding fhe commufa- 
fion befween fhem using (6.15). To consfrucf fhese and higher conserved quantifies, one may sfarf by 
simply replacing fhe field u of fhe classical conserved quanfifies by T, lake fhe normal ordered prod- 
ucls, and finally add correction ferms such fhaf fhe resulfing quanfum conserved quanfifies commute 
among themselves. 

The operator T{x) may be Fourier transformed: 

f'(rc) = y;L„e-’“-^. (6.22) 

n 

where the constant ^ is introduced for convenience. The commutation relations between the L„’s 
may be inferred from (6.15) and is the well known Virasoro algebra: 

[Lm, Ln\ = {m- n)Lm+n + “ m)5m-n- (6.23) 

The mutually commuting quantities then take the form: 


ficyl) 

-'i 



n=l 


c + 2^ c(5c + 22) c 
12 2880 M 


(6.24) 

(6.25) 

(6.26) 


The commuting variables, ^ quantization of the Korteweg-de Vries Harmiltonians H 2 k+i 

(up to unimportant numerical pre-factors) [11-13,36], in the sense that, in order to construct them, 
one may use the classical conserved quantities, which can be shown to Poisson commute given (5.23), 
and add to them quantum corrections, such as they would commute under the quantum commutation 
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relations (6.15). To fill in the pieture of how the elassieal algebro-geometrieal solutions, with their 
separated variables, 7 ^, A*, are quantized, one must show that the quantum variables 7 *, A* separate 
the problem of simultaneously diagonalizing the Hamiltonians ease, and the 

mutual eigenfunetions of have a faetorized form: 


N 

V’(7i,72,---,7Ar) = WQjilj): (6-27) 

i=i 

^{cuD 

while the eigenvalues of the operators the separate wave-funetions Qj may be found using 

the algebraie Bethe ansatz. We shall, however, not delve into this subjeet or provide any of the 
teehnieal details, whieh is a subjeet of a large body of work, whieh for the speeialized problem of 
quantum Korteweg-de Vries was first treated in Refrs. [11-13]. The message of this seetion being 
primarily that the algebro-geometrieal solutions may be viewed as the semielassieal analogues of the 
states produeed by the algebraie Bethe ansatz. The limited set of teehnieal details provided here serve 
merely to show on what basis sueh a eonneetion ean be made. 


6.3 Form factors 


We now introduee the objeet that will be eentral to the proposition of how to eonneet the elassieal 
problem of Laplaeian growth to eonformal field theory. This objeet is known as the ’form faetor’. To 
introduee it we shall eonsider the transformation properties of T. 

The eommutation relation (6.15) eontinue to hold if we replaee T{x) by a~‘^T{ax) and T{y) by 
a~‘^T{ay) . In faet for any analytie funetion z(x), the eommutation relations are invariant under the 
substitution of T{x) with T{z) with 


f{z) 


z'{x) 


-2 


f{xiz)) 


c 

12 


{z\x] , 


where {z\x} is the Sehwarzian derivative 


{z\x} 



(6.28) 


(6.29) 


This faet allows one to define an operator f{z) on the Riemann sphere (namely on C U { 00 }) by 
applying the exponential map z(x) = e*^. We obtain: 

f{z) = z-2 (f(log(z)) + . (6.30) 

Thus, for example, we obtain: 

ficyl) ^ ^ i 

^(cvl) 

The eommon eigenstates of primary fields of eonformal field theory and their deseen- 

dants. 

There is another way to obtain eonserved quantities on the Riemann sphere, whieh yields different 
eommuting quantities. Note that the proof that eommuting quantities relies only on the form 

of the eommutation relations and on the faet that the integration eontour, whieh defines the eonserved 
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Figure 5: The contours of integration of the conserved quantities and The contours of 
integration on the cylinder or the complex plane are indicated by arrows. 

(cvl) '' ^ 

quantity, is closed. This means that if we take the quantities in replace dx by dz, ThyT 

and the integration contour around the cylinder by a closed contour in the complex plane surrounding 

the origin, we will get a mutually commuting set of quantities. We denote these conserved quantities 

"fC') 

by l 2 k+i- This procedure is demonstrated in Fig. 5. 

As mentioned above, the mutual eigenvalues of the commuting operators primary 

fields and descendants of conformal field theory. To find possible spectra and the properties of these 
eigenstates, one uses the representation theory of the Virasoro algebra, Eq. (6.23). We shall find 
opportunity in the next section to quote some results from this representation theory, but we will not 
review it here. 

The diagonalization of follows a different route. The reason for the different approach 

''{cut) ^{cvD 

in comparison to F^fc+i 1^ ll*® f^^l while ^^"0 scale invariant, the operator I^k+i scale as 

j^-{i+ 2 k)^ where L is sum large scale cut-off. For example, if we think of C as the cylinder of radius 
L in the limit L —)■ cx), we expect that eigenvalue of to scale to zero in the limit L —)■ cx) as L~^, 
while becomes a universal number independent of L in the limit. The latter behavior (relating 
to is more tractable using representation theory than the former behavior (relating to To 

diagonalize one rather puts the system on a cylinder of radius L (to provide a cutoff) and solves 
the problem using the Bethe ansatz [37-39]. Indeed, treating the complex plane as a cylinder with 
the radius sent to infinity, the method of the quantum separation of variables discussed in section 6.1, 
produces the set of states, which simultaneously diagonalizes the quantum Hamiltonians ^. 

Regardless of how the two bases, which diagonalize the two commuting sets T^fc+i’ 

found, one may consider the expansion of a basis vector belonging to one set, in terms of the basis 
vectors belonging to the other. These objects, called form factors, are the main object we wish to 
introduce in this section. They hold special importance for the message this paper wishes to convey, 
as we shall argue in the next section that the modulus square of such form factors are an appealing 
candidate to describe the statistics of the Hele-Shaw interface and in particular may hold the key to 
analytically obtaining the fractal properties of the interface, including the fractal dimension. 

To be able to more precisely introduce the form factors, let us first note that usually 

diagonalized in terms operators, rather than directly by states, the two points of view being equivalent 
due to the operator-state correspondence of conformal field theory. One finds operators, ^( 2 :), such 
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that: 


®p‘^(0)]=A2fc+i$(0). (6.32) 

In addition, a vacuum, |0) is defined as a special highest weight state. Namely, this state satisfies: 

Ln\0) = 0, (6.33) 

forn > — 1. An eigenstate of = Lq — ^ with eigenvalue Ai = Ai — ^,can thus be written as: 

6(0)|0). (6.34) 

The 2 : dependence of ^( 2 ) is defined by: 

[iP,^{z)]=dM^), (6-35) 

which may be more familiar when one realizes ' = L_i. In Addition, it can easily be show, based 
on Eqs. (6.22), (6.32) and the Baker-Campbell-Hausdorff lema, that the following relation holds: 

(6.36) 

As for the eigenstates of these are obtained by quantizing the Korteweg-de Vries problem, 

as described in section 6.1. In that section it was shown that the separated variables of multi-phase 
solutions of classieal Korteweg-de Vries become separated quantum operators. The eigenstates so 
obtained are then in one to one correspondence with multi-phase solutions of Korteweg-de Vries, 
which may be labelled by the hyper-elliptic surface A. We thus denote such eigenstates by |^). In 
particular: 

iiZM) = H2k+i{A)\A) (6.37) 

We may now introduce the form factor. This is given by: 

(^|4>(2)|0). (6.38) 


The form factors have had an important role in the study of integrable field theories, and a full un¬ 
derstanding of them, would allow also for a full understanding of the integrable field theory beyond 
the spectrum. Important advances in this direction have been made by F. Smirnov (see [40] and refer¬ 
ences therein) and others (r.g., Refrs. [41,42]) . How these objects stands in relations to the quantum 
Korteweg-de Vries and its semiclassical limit are described in Refrs. [43-46]. 

The importance of the form factors lies in the fact that the set of states, |^), is, by definition, a 
complete basis and as such one can write a resolution of the identity as: 


1 = 


\A){A\dA. 


(6.39) 


The form factors then appear naturally when one inserts the resolution of the identity into correlators 
of the fields $. We must note, however, that in (6.39), the measure dA must be specified. In addition, 
it is not clear that the set of states described by the quasi periodic solutions, to which A serves as a 
label, are complete. Thus, it should be assumed that in (6.39), the label A may have to be generalized 
to go beyond labelling hyper-elliptic Riemann surface. We shall, however, in the sequel, be interested 
mainly in the semi-classical interpretation of Eqs. (6.39,6.38) and in a context where such distinctions 
are less crucial. 
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7 Quantum-Classical Distribution 

We shall now discuss some properties of a certain sesquilinear of the form factor in Eq. (6.38). These 
properties will eventually clarify the reason for proposing a possible connection with the probability 
distribution P{A). To do so, we need first to consider a semiclassical version of the form factor. 
We imagine taking a wave packet of quantum states such that the packet saturates the Heisenberg 
uncertainty bound and denote the packet by \ A)^, and define fhe objecf: 

p{A) =^0|4>(0)|^^^^|d(0)|0^ . (7.1) 

Here we allow fo fake for fhe sfafe -<0|^(0) some semiclassical approximation, and fhus denofe if as 
-<0|6(0) rafher fhan as (0|l>(0). 

7.1 Time-Translation and Scale Invariance 

We would like to show that p{A) is time translation invariant, just as P{A), Eq. (1.3). Eirst, from 
(6.37) we have for any sufficiently well behaved generic (namely, one that cannot be written as a linear 
combination of operator, O: 


lim e 
5-5>0 


^ ^-tAtnpA) 


>1.4). 


(7.2) 


Now semiclassically, and by using the Whitham method, which states that any perturbation will lead, 
in zeroth order, to the Whitham evolution, we may write: 


lim e 
< 5 - 5>0 


-tAt(iP+5d] 


1.4) ~ |.4^V 


Assuming that the object: 


(0|$(0)e 


is smooth as —)■ 0, we obtain, by employing (7.2): 

which when compared with Eq. (7.3) gives: 

p{A)=p{A^^) 


lim -<0|4>(0)e 
( 5—>-0 


'<l(0)|0A=p(.4), 


(7.3) 


(7.4) 


(7.5) 


(7.6) 


To show that p also has the requisite scale invariance property, we must first consider the correla¬ 
tion function: 


(0|l>(0)<l(z)|0) = 


c 

^2AT 


(7.7) 


where 


C = (0|$(0)$(1)|0). 


(7.8) 
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this equality being easily obtained by making use of (6.36). On the other hand by employing the 
resolution of identity, Eq. (6.39), one obtains: 

(0|^(0)$(z)|0) = J (0|d(0)|^)(^|d(z)|0)dA (7.9) 

This equation may be written through p{A): 

(0|d(0)l>(z)|0) = J (7.10) 

whieh when eombined with (7.7) yields: 

p{A)dA^ d{Hi{A)f^A (7.11) 

where Hi{A) is the dispersionless Korteweg-de Vries first Hamiltonian, Eq. (5.12), for the Riemann 
surfaee A. If we interpret the algebraie Riemann by means of the Hele-Shaw interfaee, we have the 
relation Hi{A) ~ VR sueh that we may write: 

p{R)dRr^dR^K (7.12) 

We also make a note on the measure, dA. We should demand that this measure is time trans¬ 
lation invariant. We believe that this is indeed the ease, as the measure of integration for the semi- 
elassieal wave-funetion is often assumed to be the Eiouville measure [43,44,46]. This measure is 
time-translation invariant, by a well-known theorem. The faet that it is time translation invariant with 
respeet to the Whitham evolution, may abe surmised by noting that Whitham evolution is the zeroth 
order limit in any well-behaved perturbation. In the zeroth order the time translation invarianee of the 
metrie must be reeovered, and so the measure must be invariant with respeet to the Whitham evolution 
as well. 

7.2 Some Results from Virasoro Representation Theory 

Both the proposed time translation invarianee and the seale invarianee of p{A), Eqs. (7.6) and (7.12), 
eompare well, with the properties expeeted from the Hele-Shaw probability funetion Eqs. (1.3) and 

''(cvl) 

(1.5), respeetively, if we identify the ^ eigenvalue assoeiated with namely, Ai, with the fraetal 
dimension, D. It turns out, that the representation theory for the Virasoro algebra (the algebra given 
in Eq. (6.23)), gives a diserete speetrum (again under some assumptions) for the possible values of 
Ai, and henee for the fraetal dimension D. 

At this point, it is perhaps advisable to reiterate that a more satisfaetory eonneetion between 
Eaplaeian growth and the representation theory of eonformal field theory, would start from the two 
dimensional Toda lattiee hierarehy rather than the Korteweg-de Vries hierarehy. As it stands, the 
relation between the Korteweg-de Vries limit of Eaplaeian growth and eonformal field theory may not 
be satisfaetory in order to make even eonjeetural statements. Thus, the purpose of the eurrent exereise 
linking the fraetal dimension, D, to the eigenvalue, Ai, in eonformal field theory, is give the gist of 
how the eomputation may work within a framework of more eomprehensive theory. 

After these words of eaution, let us remind the reader that the representation theory of the Virasoro 
algebra, Eq. (6.23), admits representations whieh depend on c, whieh appears in (6.23) on the right 
hand side. Among those representation there are distinguished values of c, at whieh the espeeially 
simple representations exist, the field theories at those values of c are ealled minimal models. These 
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representations are simple in the fact that finite number of fields satisfying condition (6.32) can be 
chosen fhaf generafe fhe entire represenfafion. These are called primary fields and are labelled by fwo 
number r, s. We write for such a field. We concenfrafe here on such represenfafions which are 
unifary. For unifary represenfafions, fhe possible values of c are given by: 


c = 1 — 


6 

m{m + 1) ’ 


(7.13) 


where m is an infeger larger or equal fo 3, m > 3. Wifhin fhese fheories and for fhe operators 
fhe labels r and s can lake fhe values l<r<m, l<s<r. The value of Ai for which we 
denote by is given by: 


A 


(r,s) 

1 


[(m + l)r — ms]^ — 1 
4m(m + 1) 


(7.14) 


Ofher solutions of (6.32) have a dimension Ai which may be larger by an integer, Ai = + n, 

for some re > 0. 

We have focused on fhe unifary models, since fhey seem fo have fhe kind of properfies which are 
needed fo have a consisfenl descripfion of fhe modulus square of fhe form facfors as fhe Laplacian 
growlh probabilily dislribulion. Bui of course. Ibis slalemenl musl be laken wilh a grain of sail. Al 
any rale, one can focus on fhese models and see whelher fhe fracfal dimension D which is observed 
in Laplacian growlh [1,2], D = 1.71 ± 0.3, is reproduced in fhis approach as D = + re, for 


D=2 
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Figure 6: The speclrum of dimensions, D, for fhe firsl non-lrivial unifary minimal models. Only 
dimensions belween 1 and 2 are shown. Emply circles represenf firs! descendenl of primary fields 
(namely, fields wilh dimensions D = + 1), higher descendanfs do nof have fracfal dimensions 

belween 1 and 2, filled circles represenf primary fields (namely, fields wilh dimensions D = A^''’®^). 
The minimal unifary model wilh m = Q, corresponding to c = 6/7, conlains a field wilh dimension 
12/7. This field is fhe (5,3) field in fhe Kac classification. The firsl descendent of fhe Kac (5,4) 
field has fhe same dimension. Only a few fields, and fewer primaries, are above fhe line 77 = 3/2 
corresponding to Kesfen’s bound [47], D > 3/2. 
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some admissible values n, r and s. The observed fraetal dimension have long been reeognized as 
being equal io D = 12/7, within the aeeuraey of the observation. As shown in Fig. 6, sueh a fraetal 
dimension appears at m = 6 for r = 5, s = 3 and n = 0 and for r = 5, s = 4 and n = 1. 

Many other fraetal dimension, D, ean be found for, say, the first four unitary models, 3 < m < 6. 
However, given that the fraetal dimension of eluster embedded in two dimensions must be between 1 
and 2, and moreover given Kesten’s work who showed that 3/2 < D < 2, only a few possible fraetal 
dimension appear for the first four unitary minimal models, and D = 12/7 may be found among their 
numbers. We should note that Kesten’s bound was found for the elosely related proeess of diffusion 
limited aggregation [48], believed to be in the universality elass of Laplaeian growth [1]. We should 
also note that the number of fraetal dimensions respeeting the known bounds, inereases rapidly as m 
inereases, sueh that without further insight into how to determine the eorreet eonformal field theory, 
labelled by c, from whieh to take the fraetal dimension, the approaeh will not yield unambiguous 
results, this of eourse assuming the validity of the approaeh in the first plaee. 

8 Conclusion 

We have presented a possible eonneetion between Laplaeian growth and eonformal field theory, whieh 
relies on the elassieal and quantum integrability of the problems, respeetively, after reviewing some 
neeessary material regarding the elassieal integrability of the Laplaeian growth problem and the quan¬ 
tum inverse seattering method applied to eonformal field theory and the quantum Korteweg-de Vries 
equation. It is perhaps advisable to stress again the tentative nature of the proposed eonneetion. At 
the same time, it may be reasonably said that the proposed eonneetion seems to be offer an inter¬ 
esting paradigm for how non-equilibrium integrable systems may be treated analytieally. Namely, 
that a sort of speetral expansion of field theory operators may provide the probability density for a 
non-equilibrium system. The prospeet of sueh an approaeh is appealing beeause it allows to deviate 
from the more eommon analytie deviee of proposing Gibbs or Gibbs-like measures for the solution of 
non-equilibrium problems (for example, in problems of large deviations in whieh probability densities 
ean be understood and eomputed as the exponent of some non-equilibrium aetion), while allowing for 
progress within an analytie approaeh. 

A The Baker-Akhiezer Function: Existence and Explicit Eorm 

We give here an expression Baker-Akhiezer funetions defined on some Sehottky double. To do so it is 
neeessary to deseribe the Sehottky double by the ‘Jaeobi variety’. This deseription is very familiar in 
the ease of a genus 1 surfaee. In that ease There are two interfaees and two eyeles a, h. (see Fig. (3)). 
a genus 1 surfaee has the topology of a torus and a torus is well known to be representable in terms of 
a reetangle with opposite sides identified. This latter is the Jaeobi variety for a genus 1 surfaee. The 
extension of this two higher genus surfaees, when there more than two interfaees, is deseribed below. 
The eonstruetion is then used in order to write an explieit formula for the Baker-Akhiezer funetion on 
the Sehottky double. 

A genus g surfaee admits g independent holomorphie differentials. A natural basis for sueh dif¬ 
ferentials is given by Q required to be holomorphie one-forms satisfying: 



41 



Let us consider the Abel map [22,49,50] u{z) : V —)■ C®/where 


Bkj — 



(A.2) 


The Abel map is defined by: 

Uj{z) = f (j mod V, V £ (A.3) 

J ZO 

where zq is some point on the Schottky double satisfying t(z) = z. Namely, u(z) is a vector whose 
j element is (j, whereby the freedom in drawing the contour of integration to surround m times 
around cycle and n times around cycle bi leads to an additive ambiguity in the value of u of 
the form mik + nBei. To remove this ambiguity u{z) is considered to take values in the Jacobian 
CS/(Z5 + B^9). For g > 1 the image of the Abel map is a two dimensional surface within the 
Jacobian. 

The matrix B turns out to be symmetric and, because of the existence of an anti-holomorphic 
involution, B is purely imaginary. Furthermore, the matrix B is positive definite. This allows to 
define fhe Riemann fhefa as follows: 


0{u) = exp 2m 


neZ9 


n ■ u + 


fi^Bn 


(A.4) 


The Riemann Thefa funclion is quasi-periodic on fhe Jacobian: 

e{u + ek) = e{u), e{u + Bej) = (-^ 5 ^ 


The funclion ilself is, Ihus, not single valued on fhe Jacobian. Neverlheless, if may be used lo conslrucl 
all single valued holomorphic funclions on fhe Riemann surface by faking proper combinalions of if. 
For example, fhe logarifhmic derivative of 9 is already single valued. 

The Baker-Akhiezer funclion can be wriflen in lerms of fhe Riemann fhefa funclion, as well. 
However, fo wrife fhe explicil expression, we shall firsl need, for each k lo define Iwo gr-dimensional 
vectors as follows : 



(A. 6 ) 


In addilion, lei us define a sel of poinfs 7 ^*^ as a sel of g poinls on fhe Schollky double, such lhal 
7 ^*“ = t( 7 “), and such lhal Ihere exisl a differenlial dfl wilh firsl order poles al oo± wilh residues 
±1, respeclively, and wilh zeros al Ihe poinls and 7 “. We lei = uiyyf), respectively. 

We denote from here on 


Hq = H/2, to = t = n. 


(A.7) 


The explicil expression for Ihe Baker Akhiezer function, r/;+, and a function V’ , which may be 
Irealed as ils adjoinl, reads as follows: 






4>± tk{B^t + 4) 


u{z) 


(A. 8 ) 
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where <h^(rZ;|tt) is a funetion of two eomplex ^r-dimensional veetors, w and u. The funetion takes 
the following explieit form: 

u±w-J2iufi)+K 


(A.9) 


<f>^(t(;|u) = A^{w)e' 


±1 


^±2mu-B ^Im(i«) 


9 


where K is the veetor of Riemann eonstants, whieh is given by: 


Ki = 



(A. 10) 


Expression (A. 8) is single valued on the Riemann surfaee even though, as noted before, the 
Hamiltonians are multi-valued funetion. In faet, the explieit form for the Baker-Akhiezer funetion 
is eonstrueted sueh as the quasi-periodieity of the Riemann theta funetion exaetly eaneels the quasi 
periodieity of the exponent of the Hamiltonians as they appear in (A.8). 

The amplitudes, A^, are, for now, arbitrary. The faet that the Baker-Akhiezer funetion defined in 
Eq. (A.8) indeed solves the speetral problem, Eq. (4.18), will be shown in appendix B based on the 
following analytie properties that this funetions obeys: 


• As z —)■ oo'i', the Baker-Akhiezer funetions have following expansion: 


V’n = 






-yk 

i >0 ^ 


(A. 11) 


while as z —)• ooj,, the Baker Akhiezer funetions have the following expansion 

•0^ = ^fe>o ^ (A. 12) 

i>o ^ 

The parameters are funetions of the times, tk,tk- 

• The funetions have poles at the points jf, respeetively. This faet is based the Riemann 

theta funetion partieipating in the definition of has zeros at those points [22,49,50]. 

• The amplitudes may be ehosen sueh that the Baker-Akhiezer funetions are normalized by 

I -E 

requiring to be real and positive along with the following eondition 

= “ryf = TTq- (A.13) 

■^On ^On 


These analytie properties ean be aseertained from the explieit form of the Baker-Akhiezer fune¬ 
tion, Eq. (A.8). 

It is important to note that the Riemann theta funetion has g zeros, and thus the Baker-Akhiezer 
funetion vanishes at g points on the Riemann surfaee. As diseussed in Ref. [33] and reviewed in 
seetion 5.3, the z eoordinate of the points on the Riemann at whieh the Baker-Akhiezer funetion 
vanishes, ean be eonsidered as set of g dynamieal variables. One may take the Bloeh multiplier 
assoeiated with the Baker-Akhiezer funetion on the opposite sheet of the Sehottky double as another 
set of g dynamieal variables. Together we have 2g dynamieal variables, whieh, as diseussed in seetion 
5.3 in the eontext of the Korteweg-de Vries equation, are separated variables . 


43 



B The Baker-Akhiezer Function: Definition and Properties 

Assume a Schottky double, with anti-holomorphie involution, r. We assume that as z —)■ oo, the 
funetion r(z) may be expanded as r(z) = + 0(z'^“^), for some positive integer, J, whieh one 

should assume to be large. 

We have defined the Baker-Akhiezer funetion as a funetion satisfying (4.18). We need a somewhat 
more explieit and detailed definition of the speetral problem. We write this definition here as follows: 


TTl 

'4’m^rnn — Z'lp^ 
in 

(B.l) 


, = T(z)V’n 

(B.2) 


(0 


m 

m 

(B.3) 


^1'4’n ~ ~ 'y ^ Lnlnlp^, 

(B.4) 


m m 


Equations (B.3), (B.4) should be understood as written on the front side, and may be analytieally 
eontinued to the baek side. Thus for example z in (B.3) is to be replaeed by r(z) on the baek side. 

Given that the analytie eonditions the Baker-Akhiezer funetions obeys, detailed in appendix A, 
one may prove that these funetions obey a list of properties enumerated below. The eentral item in 
the list is the faet that the Baker-Akhiezer funetion solves the Lax speetral problem. The properties 
enumerated as 1- 4 below may be treated as lemmas helping to prove this eentral faet, whieh itself is 
enumerated as property 5 below. 


1. Uniqueness: 

A funetion satisfying the same eonditions as V’n’ Eqs. (A.11-A.13), is equal to A 
funetion satisfying the same eonditions as exeept the normalization eondition, Eq. 
(A. 13), is proportional to 

2. Orthonormality: 

where dfl is a differential with first order poles at oo± with residues ±1, respeetively, and with 
zeros at the points and 7 “. The differential is normalized sueh that its integral over any 
cyele is purely imaginary. 

3. Completeness: 

Consider funetions, /^, on the Sehottky double meromorphie everywhere exeept at 
where they have the expansions: 


X eS^>o z ^ oot 


(B. 6 ) 
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with poles at 'yf, respectively. The completeness property states that such functions are ex¬ 
pandable in terms of Baker-Akhiezer functions. Namely, there exist coefficients such that: 

n-t 

= E . <B.7) 

3='n.t 

where = mijr, = mij!", = mj, respectively. The expansion (B.7) holds if 

> n'l, otherwise /"'■ = 0. Similarly, /“ = 0 unless > n||". 

4. Reflection: 

V'n (r(2:)) = V’n (^)> (B-8) 

— dQ{T{z)) = dQ{z) (B.9) 

5. Spectral problem: 

The Baker-Akhiezer functions, , are wave-functions associated with the Lax spectral prob¬ 
lem. Namely, Eqs. (B.1-B.4) hold, with (4.4). 

We now give brief, semi-rigorous ‘proofs’ of these properties 

Proof of 1. 

If two such functions ipn and ipn exist, then their ratio ^ is a meromorphic function which has 

'^n 

zeros at the poles of ipn and the zeros of 'ipn- The Riemann-Roch theorem [22,49,50] states that, 
generically, a function cannot be found with arbitrarily chosen positions for the zeros. Thus, the zeros 
and poles of pj and coincide, namely ^ is entire on the compact Riemann surface, and thus ^ = c, 

'^n '^n 

where c is some constant independent on z. The normalization condition (A. 13) fixes fhis consfanf fo 
be 1 . 

Proof of 2: 

Consider fhe infegrand ippp{{t'f,}, 2 ;)y)“({ffc}, z)dQ has no singularities excepf af where if 
has fhe form z^^~^^~^^Q'"^dz and —z^"^~^'>~^^Q’'^dz, respectively. If n > m we may deform fhe 
contour of infegrafion fo surround oo^ and obfain fhe fhe integral is zero and if n < m we may deform 
fhe contour of infegrafion fo surround oo_ wifh fhe same conclusion. If n = m, we may deform fhe 
contour of integration to surround oo+ to obfain fhaf fhe integral is given by 1, by using Eq. (A. 13). 
If for n = m we deform fhe confour of infegrafion around oo_ insfead of around oo+ and apply fhe 
condifion fhaf fhe resulf musf nol depend on fhe choice of how fo deform he contour, one obfains 
'^(7,]'*' = ^4 • The latter relation nafurally joins fhe normalizafion condifion Eq. (A. 13), and so we 

record here fhe full normalizafion condifion, which includes fhis relafion, for fulure reference: 

= ^ = ^ = • (B.IO) 

?0n ?0n 

Proof of 3: 

Eel us focus on only. The proof wifh regards fo f~ goes along fhe same lines, so we shall 
omit it. 

Eet us assume first that > mj*". Given the expansion of around oo'i-, Eq. (A. 11), one sees 
that one may cancel the coefficients in (B. 6 ), one by one for j = mp ..., -|- 1 by adding 

to /+ a linear combination of 'ipp and this without affecting the leading order behavior around 004 ,. 
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_l_ ml 

More formally, there exists a set of numbers {aj } sueh that the following expansion for the 


funetion ■ip = f'^ — V'/ holds around 


(ct + 0(l/z)) 


z-'^i e-T.hz’' 


OO'l' 


(c^ + 0{l/z)) oo^ 


(B.ll) 


for some and C|. Sinee %p satisfies all the analytie requirements defining fhe Baker-Akhiezer fune- 
fion V’^+j exeepf perhaps fhe normalization eondifion, if follows from uniqueness (properly 1 above), 

lhal Ip musl be proporlional lo ip'^ +. Lei us denole fhe eonslanl of proporlionalify by a~^ ^ and wrife 


'^ = 




j=mp+l 



from whieh (B.7) follows immediafely for /+ and wilh nf = m'jp and = m^. 

If < m~lp, Ihen if follows lhal /+ obeys all fhe analytie eondilions required of ip~^+, exeepf 

fhe normalization eondifion. Indeed, fhe expansion of /+ around 004 ^ is given by fhe form on fhe righl 
hand side of Eq. (A. 12), while fhe expansion around oO'|-, is given by fhe righl hand side of Eq. (A. 11) 
wilh = 0, and fhe exislenee of a normalized Baker-Akhiezer funelion. Thus, due lo uniqueness, 
/+ = a'ip ^+, for some a. The fael lhal a = 0 follows from = 0, showing lhal if , Ihen 

/+ = 0 . 

Proof of 4: 

If one eompares Ihe asymplolie expansion of 'ip~, given in (A.11),(A.12) wilh Ihe expansion of 
'iP+{t{z)) one sees lhal Ihey maleh, based on (B.IO), and Ihe fael lhal p{z) maps lo 00 

respeelively. The loealion of Ihe poles maleh as well. The uniqueness of Ihe Baker-Akhiezer funetion, 
Ihen dielales lhal Ihe Iwo funelions musl be equal, proving Ihe firsl equably in (B. 8 ). The proof of Ihe 
seeond equably follows mueh Ihe same lines. 

The equably, (B.9), may be proved by eomparing Ihe analytie properties of Ihe left and righl hand 
sides. The poles and residues are Ihe same. The a and b eyele integrals are purely imaginary on bolh 
righl and left hand side. Two differentials having Ihe same singularities normalized lo have imaginary 
periods over any eyele are equal, sinee Iheir differenee is a holomorphie differential wilh imaginary 
periods. The laller is neeessarily zero, sinee Ihe requiremenl of zero imaginary period represenls 2g 
real linear eonslrainls on Ihe 2g dimensional spaee of holomorphie differentials (viewed as a veelor 
spaee over R). 

Proof of 5: 

The analytie funetion on Ihe Sehollky double, whieh on Ihe fronl side is given by z is given on Ihe 
baek side by t{z), whieh as 2 : —)■ 004 , has Ihe expansion t{z) = z'^ + 0{z'^~^), Ihus Ihe fronl-side 
funetion zippp, when analylieaby eonlinued lo Ihe baek side, has Ihe expansion: 


2:V’n 


+0(1/2)) 2^00t 


(B.12) 
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By the completeness property (3 above), we obtain that there exist coefficients which are non¬ 
zero for j = n — J,..., n -|- 1 and vanish otherwise, such that: 


n+l 

j=n-J 


(B.13) 


This defines fhe mafrix L wifh fhe properfy (4.3) and proves fhe firsf equalify in (B.l). By orfhonor- 
malify (Eq. (B.5)), we have: 


Lnj = 7^^ 4) Zi^nd^- 


C 


27^^ 

Similarly, considering z'ip~, fhe complefeness properly shows: 

n+J 

= Y1 ’ 

j=n-l 

and orlhonormalily gives: 


(B.14) 


Lnj ^ T'n zT’j'd^ = Ljn, 

such lhal L = L*, and fhe second equably in (B.l) is oblained. 

The same considerations allow one lo wrile: 

n+J 

j=n-l 

wifh 


(B.15) 


(B.16) 


Lnj = T{z)'ll7Zdn 


(B.17) 


(B.18) 


We nole lhal fhe analytic confinualion of fhe fronl-side function t{z) is equal lo z on C. Thus faking 
Ihe complex conjugate of Ibis equation and, in addition, using Ihe reflection properly (4), we oblain: 


Lnj — (, (t) V’l Ztp^ dn — Ljni 


2m r 


(B.19) 


such lhal L = L^. Proving fhe firsf equation in (B.2). The proof of fhe second equation follows fhe 
same lines and we omil if. 

To prove fhe firsf equalion in (B.3) we examine the expansion of O/V’n around oo^ 


= 


^n+l^Y.tkz'‘ 
z-ne-^L2yQy+^ 

Thus we have from completeness: 




I—1 ^jn I 

j=o + 


•jn I 

- - 1 -^ 


Z —)• OO'l' 

z —)■ 004 , 


(B.20) 


-diCi 


+4 


SOn 


(B.21) 


j=n+l 
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which entails: 


n+l 


di'ipt = Y1 ’ 

j=n 


(B.22) 


with 



dii^t 

^+)4' 

^On 

0 


n + I > j > n 
j = n 
otherwise 


(B.23) 


To obtain the expression (4.4), for compare (B.20) with the expansion of to obtain: 

n—1 


dlTpn - - V’n 


dl^On ^ 

t+,t 

%n 


Y1 

j=n-J 


(B.24) 


for appropriately chosen aj. Applying the orthonormality property, Eq. (B.5), to this equation, while 
taking into account Eq. (B.23), one obtains: 


L^ = < 








n + I > j > n 
j = n 
otherwise 


(B.25) 


Comparing the element L^n in (B.23) and (B.25) and taking into account 
follows from Eq. (B.IO), one obtains 


di^ 


+,t 

On 




+,t 

On 


'^On 


which 


^nn = ^ X ^ ^ i’n • (B-26) 

But (B.l) one has ^ = {Lf-)nj, which together with (B.25) and (B.26) proves (4.4) . 

The remaining equalities in (B.3) and (B.4) are proved by repeating the same methods used above. 
The only new ingredient is that one differentiates the orthonormality condition (B.5) to obtain: 


^ ^ ^ (B-27) 

which turns out to be useful in deriving the second equalities in (B.3) and (B.4) from the first equalities, 
respectively. We thus omit these further proofs. 


C Algebro-Geometrical Whitham Averaging Method 

We now describe the algebro-geometrical version [28] of the Whitham averaging method [24,25]. 
The method requires thinking of the variation of time as coming in two sorts, long and short. This can 
be formalized by the two-timing approach, which we overview in the next subsection, this followed 
first by a treatment of the averaging procedure and finally by a derivafion of fhe Whifham equations. 
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C.l A two-timing approach 

Due to (B.14) we have the following expression for in Eq. (4.9): 

4‘V) = ^ (C.l) 

where in the last equality we have made use of (A.8) and we have ineluded explicitly the to = n 
dependence of on the index, writing The expression on the right hand side of this equation 
depends only on which themselves depend on the times {f^} only through the variable, ru({ffc}), 
given by 


^({ 4 }) = E “ 1 “ tk^k-f 

k 


(C-2) 


where 


ujk — Bujk + tSk- (C.3) 

It may also be ascertained that a^^^(n) is bounded and that the functions and thus respect 
the periodicity of the Jacobian, w ^ w + ei, w ^ w + Bsi, as noted in Appendix A. We shall call 
bounded functions which depends only on time through w{{tk}) and which respect the periodicity of 
the Jacobian, ‘multi-periodic’. 

Due to multi-periodicity it is convenient to define [riJ] as being equal to w modulo the periods of 
the Jacobian, such that [n;] lies in the primitive cell of the Jacobian lattice, l^ej + ’EBij. A multi- 
periodic function is then, in fact, a function of [n;]. Now consider some initial times and then 

consider what happens to [tu] as tk , for given k, continuously increases, by an amount which we shall 
denote as Atk- Namely, we consider [ru({fj})] for tj = fj + SjkAtk- It is easy to see that,[rZ;] will 
cover ergodically the primitive cell of the Jacobian lattice, unless iuk happens to be commensurate with 
one of the periods of the that lattice (which is not to be expected generically). Thus a multi-periodic 
function is a function from the Jacobian to C (or IR), the argument of which generically covering 
ergodically the Jacobian with the increase of any of the times. 

The ergodic property will be important in what follows, as it will allow to replace time averages 
by averages over the Jacobian. Multi-periodicity and thus ergodicity hold for the matrix elements of 
L (since they are expressible through the multi-periodic ak{n)) and for any of the operator L^^\due to 
the way the matrix elements of are expressible through the matrix elements of L, Eq. (4.4). We 
can thus write: 


L = L{w{{tk]),p,I), (C.4) 

where I denotes the Riemann surface on which the function the Baker-Akhiezer function is defined, 
from which L is derived. The inclusion of the dependence of L on p is made here for future conve¬ 
nience. It connotes the form (4.9), which includes this derivative operator. The convenience for the 
inclusion of I in (C.4) will also only be revealed below, where we shall want to consider different 
multi-periodic solutions. Indeed, since for every Riemann surface of the form considered above we 
have constructed a multi-periodic solution, we may consider different solutions by considering differ¬ 
ent Riemann surfaces. The matrix L thus depends on the Riemann surface chosen. This is encoded 
in the dependence on the somewhat abstract parameters I. To make I less abstract, one may use the 
coefficient of the polynomial in two variables Q{z, S{z)) = 0 as the data contained within I, as this 
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data uniquely determines the Riemann surfaee, by specifying the complex curve underlying it. We 
shall not, however, need such an explicit representation of I. 

The Whitham theory considers finding approximate solutions to integrable nonlinear equations by 
taking slowly modulated multi-periodic solutions. To define what constitutes as ‘slow’ as opposed 
to ‘fast’, one must first determine a scale for ujk- Let fix is fhe fypical time associafed wifh Uk'- 
l^fcl r\j . Choose a small number e. A long lime scale, r, is one such fhal if is of fhe order or greater 
fhan fhe scale £~^Tk- Namely, r > £~^Tk- A shorl time scale is a scale r such that £~^Tk ^ r. Often 
one considers short time scales such that £~^Tk ^ t ^ Tk- On such time scale ergodicity takes hold, 
but time is still considered short. 

To make use of the separation of time scales formally, it is useful to utilize the two-timing method. 
In this method one replaces the time variables as tk by tk + £~^Tk, and considers only situations 
where both tk and Tk vary on the scale Tk or much larger than Tk but definitely smaller than £~^Tk- 
The variation with respect to Tk are then considered slow, and thus Tk is called the ‘slow time’, while 
tk is the fast time, corresponding to the fast variations with respect to it. A modulated solution is one 
in which the Riemann surface is allowed to vary with the slow time, Tk, while the phase w depends 
on the fast time, tk ■ 


L{tk + e-^Tk) = l{w {{tk}), P, I{Tk)) , 


(C.5) 


where 


P = 


d 


d 


dto OTq 


(C.6) 


The matrix elements of such an L are slowly modulated solutions of two dimensional Toda lattice. 

The form (C.5) conveys the fact that the multi-periodic L varies only on the long time scale, Tk- 
This form is an approximate solution to the two dimensional Toda lattice equations, which becomes 
exact in the limit e —> 0. Indeed, in this limit a variation of the argument of L by dT is accompanied 
by variation of the argument of I by edT —)■ 0, and the multi-periodic wave becomes unmodulated. 
The form (C.5) is, however, rather naive, in that it is not a good basis as a zeroth order term in an 
expansion in the small parameter e. A better ansatz is given by: 


L = L{e-^W{{Tk]),P,I{Tk)), 


(C.7) 


with 


dW _ ^ dW 
&Tk~'^^' Wk~ 


(C.8) 


This ansatz also coincides with an unmodulated solution in the limit e —)■ 0, but takes better account 
of the change of the phase of the multi-periodic wave as the wave is modulated. 

Note that the integrability of equations (C.8) imply: 


dujk _ 

dTi ~ dTk’ dfi ~ dTk’ 

and the complex conjugates of these equations. These equations called ‘conservation of waves’ con¬ 
stitute constraints on the way the Riemann surface may change. These constraints are not enough to 
determine the dynamics, but shall be useful below to obtain the modulation equations in full. 
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We shall need two more ingredients in order to proeeed. First, we shall want to inelude the first 
order eorreetion term to (C.7), whieh we shall denote by eLi. Seeond, we shall need in the sequel 
derivative of Lq with respeet to the slow times whieh appears in I only. Namely, negleeting the 
derivative of the slow times in VV. It is thus useful to introduee the times as follows 

L = Lo(e-^W{{Tk}),P,l(^{fk}))+eL^ + ..., (C.IO) 

and take Tk = Tfc at the end of the ealeulation. The operator P is still given by (C.6), as the shift 
operator in (4.9) does not depend on {Tfc}, sinee it does not hit I. 

Note that, as a result of the three-timing method, a derivative with respeet to time, dk, whieh was 
taken before the slow times were introdueed and before these slow times were split into Tk and Tk, 

becomes the expression dk + £ 7 ^ + , after the introduction of the new times. 

OJfe dTfe fk=Tk 

Finally, note also that a similar expansion holds for We write this expression for future 
reference: 

L« = 4') (C.ll) 


We shall use the boldface T; to denote both times Ti and T/. Furthermore we shall write: 


d d d 


(C.12) 


To obtain (2.51) as the semiclassical limit of the zero curvature conditions, Eq. (4.6), we start 
from the latter equations and expand them in s. We have: 


'-^+eT-LV>T+e^-LV‘>] =0. 

dti dTi dtk dTk 


(C.13) 


Substituting form and L the expansions (C.10,C.ll), and P = p + ed^, and then expanding in e, 
one obtains that the zero order terms in e vanish due to (4.6), while the first order term in e yields the 
equation: 


F = G, 


(C.14) 



dif 

dTi 

dTk 


dPf^ 

dti 

dtk 


_ s tT) t (^)l 

1-^0 ’ ^0 J 




(C.15) 

(C.16) 


The Poisson bracket on the left hand side of (C.14) is defined for any fwo funcfion of P and Tq as 
follows: 


rA/^rr,\ dAdB dBdA 

{ {p, o), {p, = 


(C.17) 


and if assumed fhaf fhe operafors acf on functions of fhe fasf variables only. Nofe also fhaf P = 
^ is equal fo leading order fo p = ^one may replace P by p in (C.14). 
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The appearance of the Poisson brackets in (C.14) is related to the expansion with respect to e of 
functions of P = p + . Indeed, given A{P, Tq) and B{P, Tq), one may expand A{P, Tq) = 

A{p, To) + + 0{e‘^), with a similar expansion for B. One then obtains: 


[A{P,To),BiP,To)] = [A{p,To),Bip,To)]+e{Aip,To),Bip,To)}+ 

d 


(C.18) 


/ 

. dB 


^ dA 

+ £ ( 

A, ^ 

— 

B, 

V 

dP 


dP 


dTn 


+ 0{ep. 


Since we assume the operators act on a function of fast times only, we can drop the third term on the 
right hand side. 


C.2 Averaging 

We define the averaging of any matrix. A, as: 


(A) = lim 


no+L 




LToo{ 2 Ly ^ ^ Jl LLV 

'' ’ n=no-Li,ieZ 


Ann n dtk- 


k=l 


(C.19) 


The scale L is to be intermediate between the fast and slow scale. Namely, Le, must be thought 
as a small number. Since this small number does not, however, scale with e, for the purpose of the 
expansion in e, one may assume L~^ ~ e. Note that in the limit L —oo, the sum over n becomes a 
trace, and thus one expects: 

([A,P])~i~e. (C.20) 

For this to be true, it is necessary that the object being varied over, [A, B], is bounded within the 
domain being averaged over. 

An important average that we shall encounter will be of the form: ([foA, P]), with to being a 
matrix with elements given by: 


{io)nm — 'kldnm- (C.21) 

We shall need an alternative representation of such an average. For this purpose note that the average 
can explicitly written as: 


Uf 

([tflA, P]) — ^ ^ ^ ^ {pAn^n+jBn+j^n j)An—j^nBn,n—j)^f^^ ; (C.22) 

j n=ni 

where rii = uq — L and nj = no + L. Here the average on the right hand side, denoted by is 

an average over the times t^, with k > 0, only. This expression may be further simplified by wrifing: 


ni+j\ 

{nAn,n+jBn+j,n)^,^^ . (C.23) 

n=nj-j n=ny 

Since j <C L, one may replace n by nj in fhe firsl sum and n by n* in fhe second sum. In addifion, if 
A and P are mulfi-periodic, all fhe ferms in fhe firsl sum are equal fo each ofher in leading order, due 
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to ergodicity, and the same is true for the seeond sum. Taking all these considerations into account, 
we may write: 

{[ioA, B]) = ^ ^ J {rifAnj^rif+jBnf+j^nf ~ '^i^ni,ni+j Brn+j^m') (C.24) 

j 

which in turn can be written as: 

Uf-l 

{[ioA,B]) = ^ ((n + l)(^[fo,^])n+i,n+i -n{A[io,B])n,n)^i.^y + 0(e). (C.25) 

n=ni 

The fact that this equation is equal to the one preceding it easily follows by returning to a representa¬ 
tion of {A[io, B])n,n as J2jj^n,n+jBn+j,n, and similarly for (A[fo, B])n+i,n+i- Ergodicity implies 
that, to a high precision, {{A[io, = which leads to: 

{[ioA,B]) = {A[io,B]) + 0(e). (C.26) 

Finally, we may use the equation [fo, B] = — which holds generally to write: 

([t„A,B|) = - (a^) + 0(E) = - (^a) + 0(E). (C.27) 

where in the last equation we used the fact that the average over the commutator of two mulit-periodic 
functions is small. This equation, which holds for quasi-periodic, A and B, will be useful in the 
following. 

C.3 Krichever’s derivation of Whitham 

The algebro-geometric version, due to Krichever [28] (see also [27]) of the Whitham averaging 
method, consists of averaging the nonlinear equations with the Baker-Akhiezer function. Thus, we 
shall consider objects such as {Ail;^. Here t denotes the transpose, such that A^|J~^ 1 p~^ is a ma¬ 
trix whose mn element is i’n- Furthermore, the Baker-Akhiezer function is assumed to 

depend on the fast times through the variables {ffc} that appear explicitly in (A.8), the dependence 
on the slow times is implicit in this explicit formula, through the meromorphic differentials, Hk, the 
frequencies and indeed the definition of the Riemann theta function, Eq. (A.4), which depends on 
B, which is itself a function of the slow times. Note that 

(^+^-‘)n,n+, = (C.28) 

and following the multi-periodicity of the matrix elements of are also multi-periodic. 

Applying this averaging procedure to (C.14), we shall show that one obtains: 

0 = (C.29) 

with F and G defined in (C.15) and (C.16), respectively. We shall show this by showing that 
= 0, and then F = G gives the second equality. Indeed, consider averaging applied 
to G. One easily obtains from (B.3): 

^ (C-30) 


53 



Averaging this equation, the left hand side vanishes to leading order as an average over a full deriva¬ 
tive, and one obtains: 


0 = 


dL[^^ 

dti 


if; * 





This ean be written as: 


(C.31) 


0 = 


dL\ 


(k) 




dti 


T (^) T (0 

-^1 > ^0 






(C.32) 


As an average over a eommutator of multi-periodie funetions, the last term on the right hand side 
vanishes in the limit L —> oo (or rather is of higher order in e). If we follow the same proeedure with 
k and I interehanged, one obtains = 0 and thus = 0 

The speeifie terms in the equation = 0, when one inserts the definition of F, Eq. 

(C.15), ean be further developed to obtain the modulation equations we seek, Eq. (2.51). We start 
with the first term in the definition of F, that is the first term on the right hand side of Eq. (C.15). 
To be able to re-write the average over it with the Baker-Akhiezer funetions, one first eonsiders the 
following objeets: 


JLJL 

mWk 

- JL 
~ m 







(C.33) 


where the prime denotes taking the funetion at slow times {T'} rather than at {Tj}. We may now take 
{Tj} = {Tj} in this equation to obtain: 


d d 

WiWk 


dL^k) 

dT, 


11 )^11) * -|- 





The derivative reads explieitly: 




m>0 


dHffi dujm 



+ to 


dHo 

dTi 


dijJo ^ 

H- - ■ V 

dTi 




(C.34) 


(C.35) 


where to is defined in (C.21) and the gradient V"' denotes taking the gradient derivative of with 
respeet to its first, veetor, argument. Namely, 




w=w({t^.}) 


Inserting (C.35) into (C.32) and making use of (C.27) gives: 


d d 





* — 




dHo 


dp ydTi dTi 




-t 


(C.36) 


(C.37) 
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On the other hand taking the opposite order of the derivatives in (C.32) one obtains: 


d d 
WkWi 




A 

dtk 



dHm diOm ^w\ 
dTi dTi ' ) 



(C.38) 


Averaging over this equation, and disearding averages over full derivatives (note that one averages 
over all times and thus, e.g. {tmdk ...) = 0, by ergodieity, even for k = m.), one arrives at: 


d d 

WkWi 



dTi ^ m ) 




Equating (C.37) and (C.39) one obtains: 




ST, J 



(C.39) 


(C.40) 


giving a representation of the average over the first term in the definition of F, Eq. (C.15). In that 
respeet note that exehanging by has no eonsequenee to leading order. 

We may subtraet this equation the same equation where the roles of k and I are reversed. Making 
use of the eonservation of waves, Eq. (C.9), one obtains: 



To deal with the average of the third term in the definition of F, Eq. (C.15), one eonsiders: 


d A 

dT^ \^dti 



d 



A 

dtk 


dp 




'll; '^ip 







(C.42) 


where G is defined in (C.16). Taking {T} = {T'}, one obfains: 


I d (dL^^"> OA T-i\ 
\dti y dp dTo j 


(I-AA'l 

\dtk y dp dTo J 




(C.43) 


where fhaf faef fhaf *) = 0 for all Tq has been used. One may eompufe fhe lefl hand side 

of (C.43) by firsl applying fhe derivafives wifh respeef fo Tq and only fhen faking fhe derivative wifh 
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respect to ti or tk- Much the same methods as were applied to obtain (C.41) allow one to write: 




dHi. / 


on 


dp 




-t 


+ 


(C.44) 


dcoi / 




—t 


dTo \ dp 


dujk / dL^^^ 
dTo I dp 




-t 


We are now in a position to write down ),by combining (C.41) and (C.44) with the 

definition of F, Eq. (C.15). Taking into account the conservation of waves, Eq. (C.9), yields: 



which combined with (C.45) gives: 


0 = 




Ho- 


dHo 

dTk 




dHo 

dTi 


Hk 


(C.45) 


(C.46) 


(C.47) 


These equation must be true for all z on the Riemann surface. Examining the analytic properties of 
each term of this equation, one obtains that an over-determined set of conditions on the functions 
multiplying Hq, Hi and must be satisfied if fhe righf hand side is fo be equal fo zero, implying fhaf 
each summand is zero by ilself. Eeading fo: 


dHi dH,. 

-^-^ = 0 

dTj dTi 


(C.48) 


for any i and j (including any fwo chosen from fhe sef {/c,(,0}), which is equivalenf fo fhe firsl 
equation in (2.51). The same mefhods may be applied fo obfain fhe second equafion, namely: 


dHi dHj 
'Wj ~ ~Wi 


(C.49) 


We shall give defails of how (2.21) is obfained from (4.8) by averaging since much fhe same mefhods 
are employed here as well. 
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